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by the use of even or odd polynomials, i.e., 
sin(x) = cP (ae) 


coe(u) = Piao) 


where v. and P. are the sine and cosine polynomia! approximations. 

Table $f.A-f gives the ranges over which the various functions are 
approuimated By polynomials. The various range reductions relations used are 
Listed if the program. ‘ccc Appemtta AL) 
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TABLE II.A-1 


Ranges over which polynomial approximations are used 


for the various functions. 


Function Interval 
SIN t om/u, 
CoS zt w/u 
TAN < 4/8 
ATAN a(V2 <1) 
aS in £0. 375 
EY 0, 0.5 
uM 0.5, 1.0 
SORT 0.5, 1.0 
S{MK 21.0 
COSH 21.9 
TASH £0 .5 
AS (eH 20 . 375 
at any £0.25 
acOS polynomial not used 
acOSH polynomial not used 


taken. The operations dealing with the powers of 2 are assumed to be very 
fast; they will be machine dependent and should be tailored specifically for 
each machine. 

The hyperbolic functions and their inverses are approximated by poly- 
nomials in a narrow region around the origin and outside this region are 
computed from their relationship to EXP, LN and SQRT. 

The inverse functions arccos(x) and invcosh(x) have an anomalous 
behavior in the neighborhood of x = 1; they cannot be approximated by a poly- 


nomial in x near x = 1. Instead we have 


arccos (x)= V2y + y°/3 + Wy 3/45 + aes 
where y = 1-x, and 

inveosh(x) ~*~ V 2y ~ y°/3 + y>/45 eer 
where y = x~-l. 

The first two terms are used to approximate the functions,and the third 
term is used to control the range over which the approximation is used in such 
away as to maintain the error tolerance. Outside this range the arccos ne is 
computed from the arcsine; the invcosh is computed from LN and SQRT. 

For the arcsine, an odd power series is used in the region about the 
origin; outside this region the arcsine is computed from its relation to 


arctangent and SQRT function, i.e., 


arcsine(x) = arctan(x/ V 1-x7 ). 


II. B. Chebyshev Approximations 
The advantage of the Chebyshev polynomials is that they are the poly- 


nomials having the largest number of maxima and minima on the given interval 
and for which all maxima and minima are equal in magnitude. Thus, if the 
lowest order term neglected can be considered the error function, it can be 


seen that the error is more or less uniformly distributed over the interval, 


and can be shown to have the smallest maximum error of any polynomial approxi- 
mation of the same order. 

Since accuracy requirements are known at compile time, it is possible in 
Ada to determine doth the number of Chebyshev terms and the values of the 
coefficients to the power series at that point. 

Example: sin(#x). 

This function can be mapped onto the interval -i<xs+1. The functio can 


then be well approximated in this range by the n-term Chebyshev series 


where the Bs can be calculated from 


K 


] 
B = al T(x) sin(wx)dx/V 1-x- . 


-1 
(Only the odd numbered Bs will be nonzero.) 


Each TL Ox) is a polynomial of order k 


_ j 
T(x) = >» ae 


o<jsk 


The Cs are well known and calculated exactly. 


The n-term power series approximation for the sine then is 


sin(ax) =® »S », B,C, 4X" 


lfkEn o<j<k 


a ax? 


osjsn 


eae 


where 


> 
t 


7; » BK ks: 


jsk<n 


We see that adding more Chebyshev terms, say to improve the accuracy, changes 
all the A coefficients of the power series. Similar considerations apply to 


approximations for the other functions. 


II.C. Other Approximations 


In addition to Chebyshev derived polynomials, two other approximations 
commonly used are Pade’ like approximations and continued fraction. The 
Pade like method consists of approximating as a function a ratio of two 
polynomials P(x)/Q(x); it can be converted to an equivalent continue fraction 
and vice versa. While these methods can sometimes be superior to a simple 
linear Chebyshev polynomial approximation in the sense that the same or better 
accuracy can be obtained with fewer computer operations, they all suffer the 
same drawbacks from the point of view of trying to write machine and accuracy 
independent programs. This is that for optimum conditions, all coefficients 
or constants in the method will change if the accuracy is changed. The Pade' 
and continued fraction methods suffer other problems. 

First, it is much more difficult to find the optimum set of coefficients 


for a given accuracy, while this is relatively easy to do for the Chebyshev 


series. The Chebyshev series requires storing only a single set of constants 


(10) for each function while for the P(x)/Q(x) method, an entire set of 
constants would be needed for each accuracy interval. 
A second reason for not going with the ratio of polynomials is that 
it requires a division operation and the trend for hardware today is toward a 
relatively slow divide operation relative to addition and multiply. Present 
day hardware suggests we should avoid divisions where practical. (See IV.A.2 
for a further discussion of this point.) 
In summary, the Chebyshev method was chosen over the P(x)/Q(x), the 
ratio of two polynomials, because 
1) Variable accuracy methods are relatively easier to obtain; 
2) Fewer stored constants are required; 
3) A division operation is eliminated (or traded for some number 
of multiplications and additions). 
III. Ada Considerations 
The principal Ada techniques expected to be most useful in developing 
machine and accuracy independent functions are the package and generic con- 
cepts. The main idea is that the functions will be embedded in a package with 
a generic parameter, call it REAL, describing the type (number of digits, to 
be used in the floating point arithmetic. The type REAL or its nuaber of 
digits is then specified by the user in his main program. Then at the time the 
package and user program are integrated, the package is instantiated with the 
appropriate number of digits for REAL. 
There are a number of program features that depend on the number of 
digits in REAL (called in Ada REAL DIGITS). 
First, the number of terms in the Chebyshev series and so the power 
series is determined by requiring that the Chebyshev terms not used are all 
smaller than 10##(-REAL'DIGITS). Next, this determines the size of the arrays 


that hold the power series polynomial coefficients. Then, in some cases 


10 


eee 


(ACOS, ACOSH, and TARH) the range in which a particular approziaation i8 used 
will be determined by the number of digits in REAL; tne computation of these 
ranges isa carried out in the initialization Block of tne function pachare. 

Certain assumptions concerning the computer arithmetic have deen aade in 
writing these packages. We delieve these assumptions are, in general, consise- 
tent with the Ada view, the IEEE standards, and the MCF or Mebula choice of 
floating point arithmetic. These assumptions are as follow: 

1. There ta a gaaximm floating point nusbder (called in Ada 
FLOAT‘ LARGE). 

2.  =FLOAT*LARGE exists and is the gost negative floating point 
numbder. 

3. Every floating point number including FLOAT'SMALL, out 
excepting zero, has a reciprocal; the reciprocal of 
=PLOAT*LARGE may de zero, but need not be. 

&. Any positive number can be subtracted from FLOAT'LARGE and any 
positive number can be added to -FLOAT LARGE without causing 
an exception alara. 

5. FLOAT'LARCE can be divided by any number greater than or equal 
in magnitude to 1.0 or muitiplied by any number less than or 
equal in magnitude to 1.0 without causing an exception. 

6. Any two floating point numbers (including TFLOAT'LARGE, can be 
compared without causing an overflow or exception. 


IV. Routines 


IV.A. Introduction 


IV.A.1. Variable Accuracy 


One of the prime goals of this effort was to construct these subroutines 
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froa the corresponding Chebyshev coefficients. These Chebyomey ope!s 3: iis 
nave been precalculated and are stored a8 part of the progranm of pa t+ ne: 
containing the function. Enough teres in the Credysnrey series are sist 1: 


= 16 


be able to achieve eA accuracy (relative error) of 10 It 19 ese wet tha’ 


thie accurecy lisitation (19°, will be adequate for the wast @ejor:iy ° tte 
enticipeted Ada applications in embedded systems; estenpian tc higtes ai: -- 
recy would be straightforwerd. Table [V.A.1-1 gives a tedle of the a - +>, 
versus the numéer of teres needed for each function. There is, of comers. ‘ot 
each function @ certein amount of computation necessary to reduce tt< 2 «€ m=: 4 
renge to within the range of the poiynopial., Owt that is not refiecte = .2 i ua 
table. (Also not reflected in Table I1V.A.1-f is the fact that certai* 
funetions call other functions for same renges of their welves rather than 


weing a polynomial.) Table IV.A.1-11 lists those functions that cali others 


and in what ranges. 
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There is one exception to the above description of how the error is 
controlled and that is for the SQRT function. In the square root function a 
single Newton-Raphson or Heron's iteration is performed; i.e., 

xe (y/x # x)/2 
where xe VV . This operation squares the relative error when the error is 
gmail. For SQRT it is only the first approximation that is obtained from the 
Chebyshev polynomials, the final value is obtained from the above iteration. 
Thus the accuracy of the power series need be only as good as the square root 
of tne relative error tolerance; i1.e.,if the relative error tolerance is 107°, 
the polynomial need only be as accurate as 107°"? & 3.167> -~ the Heron 
iteration will igprove this back to 107°, 

The Chebysnev approximations normally give absolute not relative error 
ligits. When dealing with floating point arithmetic it is relative error that 
one wishes to saintain more or less uniform over some internal. For the even 
functions this is not a problem since in the interval of approximation the 
functions do not change greatiy and so the relative error ‘s very nearly equal 
to the adsolute error (the functions being approximately 1 at the origin). To 
@aintain the relative error sore or less constant for the odd functions, say 
f(x), we use @ Chedvshev approximation to f(x)/x. This gives nearly uniform 
relative error for f(x) on the interval around x = 0,and again the absolute 
error i3 approximately equal to the relative since for these odd functions 
Li@ (x0) F(e)/n 2 3. So for example, the Chebyshev coefficients for the SIN 
function are actually those for sin(x)/x, and those for TAN are in fact for 


tan(a)/x, etc., for all the odd functions. 
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IV.A.4. Error Conditions and Out of Range Alarms 


In nearly all computer arithmetic systems there is a limit to the size of 
a floating point number; there is also a smallest-in-magnitude nonzero number. 
These numbers are machine dependent but referable in Ada as F'LARGE and 
F'SMALL, where F refers to the (floating point) type. Because of these limits 
there are certain functions such as EXP, SINH, COSH, etc. for which only a 
relatively restricted range of inputs are valid. Thus for example, if we 
refer to our floating point type as REAL, input values to the EXP function 
greater than ln(REAL‘LARGE) cannot be computed since they would result in 
values greater than REAL'LARGE. There are other functions such as LN and SQRT 
that havc invalid ranges because the result would not be in the real number 
system (but would be compiex). For example, numbers less than or equa: to 
zero are invalid inputs for LN; X <0 is invalid for SQRT(X). 

The problem of what to do about invalid inputs is complicated by the fact 
that all computer arithmetic with floating point numbers involves rounding; 
floating point arithmetic is seldom exact. When rounding takes place in the 
neighborhood of the boundary between the valid and invalid region of the 
function, it is difficult to decide what is the best action to take. Thus for 


example, if one wants to evaluate: 


SQRT(A + B + C = SIN(D)) 


and the exact value of A+B+C-SIN(D) is 0 but due to rounding (or truncation in 
SIN) the computer results turns out to be -1.0 E-9, what is the proper action? 
How should the SQRT subroutine respond to small negative inputs, and what is 
small? One action would be to trip a numeric exception alarm and let (force) 


the user handle the problem. Some systems halt the computation at this point. 
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Another solution is to assume that there is a small negative region such that 
input from that region ought really be treated as input of zero and not trip 
the exception alarm for these inputs, simply return zero; however it is 
difficult to determine how large or small this boundary interval should be. 
The solution to the problem of the invalid input here is different yet 
from those mentioned above. We take the point of view that there should be no 
truly invalid input region and redefine our functions so all inputs are 
legitimate and will return some value and no exceptions are ever raised. For 


example, SQRT is defined so that 


sort(x) zeV¥% 29; 


0 X< 0 


The square root function given here always returns a value, with no exception 
raised for input values from -REAL'LARGE to REAL'LARGE. Thus, small rounding 
errors in the vicinity of zero will cause no difficulties; however, it is the 
users responsibility to trap any input values to SQRT that the user feels 
ought to be considered illegitimate. Table IV.A.4-I gives the extended 
definition of the elementary functions showing how the usually considered 
invalid regions are treated. As an example, Figure IV.A.4-I shows how the 
aresine function is extended beyond the normal input range of = 1.0. 

As mentioned above the main virtue of extending the definitions of the 
functions into their normally invalid regions is to avoid extraneous error 
signals when the argument drifts into the invalid region by a small amount of 
rounding or truncation error. This solution may not be the best resolution to 
the problem nor the solution desired by a particular user. It does however 


have the virtue of allowing the user to determine what he would like to see as 


Ej2> 


dDNV 1 , TVau O°? O°T 


——_____-. - - 


ee cfu 


(X) NISV 


O-6> ADAV 1, TVIN - 
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Table IV.A.4-I 


A list of the functions and the way in which their inputs were extended to 


the normaly invalid regions, together with the ouput values in these regions. 


SIN(x) = sin(x) all x 
COS(x) = cos(x) all x 
TAN(x) = tan(x) all x 


(returns * REAL'LARGE 


for x = 7 ®/2, *3.8/2, *6 9/2...) 


ASIN(x) = ¢arcsin (x) Ix} < 1.0 
+ %"/2 x > 1.0 
~ W/e2 x < =1.0 
ACOS(x) =< arcos(x) 1x} < 1.0 
0 x > 1.0 
~~" x < -1.0 
ATAN(x) = arctan(x) all x 
EXP(x) = 2 e* Ix} <  1n (REAL'LARGE) 
0 x < ~<ln (REAL'LARGE) 
REAL ‘LARGE x > ln (REAL'LARGE) 
UN(x) = § In(x) x > 0 
-REAL * LARGE x <0 
SQRT(x) = <¢vVx x 2 0 
0 x < 0 
SINH(x) = ¢ sinh(x) )x} < In( REAL' LARGE) 
REAL * LARGE/2 x >1n( REAL‘ LARGE) 


~REAL *LARGE/2 x < =-ln (REAL‘LARGE) 
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Dee rag 


Table IV.A.4&-] 


GC uyeart mt) 


TANH(x) = tanh(x) all x 
ASINH(x) = invsinh(x) all x 
ACOSH(x) = ¢ tnve-sh(x) x 21.0 
0 x<1.0 
ATANH(x) = ¢ invtann(x) mp< 1.0 
REAL * LARGE x > 1.0 
REAL ‘LARGE x < -1.0 
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ene valid regions and alious (requires) Bim to eet Ale om error traps for 
what tnat user considers illegitiaate inputs. 
IV.8.%. SIN ams (OS 

Input Valve Hange: -NEAL* LARGE tu MEAL ‘LARGE (in radians) 


Output Value Range: -!.0 %6 «3.9 


Tne argument, «, 18 first aapped onto the interwal -3 @/2 to 3 872. 
tf x <¢ 3/3 tne sine of cosine polynomial :s computed for tne respective 
function; if ‘3, is greater than F/8, SIN calis tne cosine polynomia)d and (as 
calls the sine polynomial with the appropriate sign end phase ad )justeent . 
There are no invalia input or output ranges for SIN and COS. However 
vecause of truncation and/or rounding errors it is possibile for results to be 
returned which are slightly larger in aagnitude than 1.0. 
IV.B.2. TAN 


Input Value Range: -REAL’LARGCE to NEAL‘ LARGE (in radians) 


Output Range: ~-REAL'LARGE to REAL LARGE 


The argument. x, is first sapped onto the interval -59/6 to 58/8. If ix 
then is less than 7/68 the tangent polynomial is calied; if ixi its in the 
range 7/8 to 39/8 the tangent polynomial for |x} - 7/8 i8 computed (as Q) 
and TAN is evaluated as (Qe!l)/(1-Q), with appropriate sign’. If txt > 3 8 /8, 
TAN is computed from the reciprocal of the tangent polynomial with appropriate 
sign and phase adjustajent. 

The only invalid input arguments are Ixl 2 9/2, asi, 3, 5, ete. For 
these angles the value = REAL* LARGE ia returned. No error halt or exceptions 
are generated. 

+ This is an application of the relation : | 
tan x = (tan(x-b) + tan b)/(I - tan(x-b)) 


with b = t ¥/4 
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IV.8.13 ATANH 
Valid Input Value Range: -!.0 to 1.0 


Jutput range: -REAL‘LARGE to REAL*LARGE 


.f tne argument, *, $8 tm tne interval -0.25 to 0.25 a polynomial 
approsimation 18 used for ATANKH. [In tne interval 0.25 <Ixl< 1.0 the relation 
a (i) @ a/v tees) /2 23 wned. 

For the inwaild input region ixi + 3.0, the value = REAL‘LARGE _> returned 


aft mo error halt of exception is raised; tnis point is discussed in Section 


Verification of the Function Paci 


Short of proving that each function given nere is correct, the best that 
7anh Se done ia to chet the sorrectnmeas °F tne vaive coeputed for a large 
mamber of arguments for each function. ence at this point we have no way of 
obtaining check valves internaliy in Ada, it appears that to check then, we 
Quest write them oul and verify each vaive by hand or by another computer 
program ina different :anguage such as FORTRAN. 

There is, however, another good method of verifying the correctness of a 
Large number of valves af a function and that is by checking certain addition 
theorems for these functions. Fer erxanple the sine and cosine should obey the 
relation sin’x ® 309" : ' for ali x. This, however, i838 not a very good check 
since the relative error of the smailer of these could be very iarge but not 
be observed. A better check are the “triple relations*. For example 

sin 3a =: 3 sin x= 4 sin?s 
can be used to verify the sine function. The triple relations have the 


advantage of requiring just one of the functions at a time and allows it to de 


3$ 


—_ ~~ — — 


checked against itself. Table IV.C~I gives a set of triple relatios for the 
functions considered here. A very good verification procedure would be to 

! check the triple relation of each function for several thousand (random) 
points uniformly distributed overa range that would ensure the exercising of 
all branches of the function code. Let us emphasize again that this has not 
been done yet due to lack of time and not having a suitable Ada compiler. Only 
the exponential function has been executed and partially verified; Figures 


V.B-f and II give triple relation error curves of EXP for a limited sample of 


arguments. 
V. Critique 


V.A. What Went Wrong 

The most serious difficulty encountered by this project was the lack of 
an adequate Ada compiler. In fact not until 30 days before the final 
termination date of the project (60 days after the original termination date, 
a three month no cost extension was requested and granted), did we obtain the 
Ada/ED compiler. Ada/ED (Version 11.4) was obtained through NTIS and is 
aonsidered a preliminary unvalidated version intended for education and 
experimental use only. And while it was extremely complete, it did have a 
variety of bugs or errors; those we uncovered are listed in Table V.A-I1. 
These er-ors, while not serious once they were understood, contributed 


significantly to our confusion. We were, in fact, learning about Ada and tne 


Jifficuities of untangling our mistakes and misapprehensions about Ada from 
~. rors ir Version 11.4 proved a real strain and lead to the consumption of 
vast amounts of computer time. (Ada/ED is an experimental version and not 


\ 
du -igned for production use -- it is very slow -- by a factor of about 10 to 
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Table [V.C-I 


"Triple Relations" for Verifying the Accuracy of the Elementary Functions. 


The 
sin 3x = 3 sin x - 3 sin?x 
cos 3x = -3 cos x + 4 cos >x 
tan 3x = (3 tan x - tan>x)/(1 - 3 tan°x) 
e* = (o*)3 
In 3x = ln 3 + ln x or ln x? = 3 in x 
V3x = V3Vx or Vid = WX )3 
sinh 3x = 3 sinh x + 4 sinh>x 
cosh 3x = -3 cosh x + 4 cosh>x 
tanh 3x = (3 tanh x + tanh>x)/(1 + 3 tanh°x) 
3 sin! x = sin” | (3x - 4x3) 
3 cos”! x = cos™!(-3x + 4x?) 
3 tan) x = tan 7) ((3x = x3)/(1-3 x°)) 
3 ginh7! x = sinh™'(3x + 4x?) 
3 cosh”! x = cosh™!(-3x + 4x?) 
3 tanh”! x = tanh”! ((3x + x )/(1 + 3x°)) 
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A list of known flaws (as of 10/26/81) in the Ada/ED Compiler Version 11.4 


Compile Time: 


1. Generic and actual parameters may not have the same name. 
2, LONG FLOAT not implemented. 


3. Won't handle exponentiation of a universal constant. 


Run Time: 
1. Won't multiply by a floating point number if its value is 0.0. 
2. INTEGER truncates instead of rounds (sometimes?). 
3. Universal constants as parameters of functions or procedures cause a 


RUN time error halt, not a compile time error. 
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10° (in execution) over what will be required in a production mode.) 

The net result of not obtaining a compiler until so late in the program 
and then to have it be somewhat flawed has been that the functions and 
packages here are not very well debugged. In fact we will have to include the 
disclaimer that “this package of functions is an unvalidated early version and 
is intended for experimental use only". We would like to be able to claim at 
this point that each of the functions has been thoroughly debugged and can be 
certified correct; but this is not at all the case. The best we can claim is 
that they have been syntax checked; except for the EXP function, none of the 
functions given here have been executed. Preliminary debugging of the EXP 
function has been accomplished and error curves obtained showing that EXP does 
meet the error tolerances requested. however, even here more checking should 
be done before it is released for incorporation into a "live" system. 

V.B. What Went Right 

The main bright note of this project is that most of the ideas about how 
to use Ada to build library packages seem to have checked out. To sumnarize 
these briefly: 

1) Packages; the packaging feature of Ada works nicely and is just what 
is needed for writing library subroutines. 

2) Generics; also work well and provide a mechanism for constructing 
variable accuracy routines where the error tolerance is effectively supplied 
by the user at the time the package is integrated into his total program. 


This has been demonstrated explicitly for the EXP function where the error 


curves show how the accuracy depends on the users specification of the number 


of digits in his REAL type variables; see Figures V.B~I and II. 
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3) Overloading; this was a useful and convenient device but was not 
really crucial to setting up the library routines. The main advantage of 
overloading was that it allowed mixed mode arithmetic (most of the arithmetic 
operators were overloaded for various combinations of INTEGER, REAL, FLOAT 
and LONG FLOAT). 

4) Function and procedures; these constructs also worked well and of 
course were indispensable in setting up an elementary math function library. 

5) Dynamic Arrays; worked well and were especially useful for library 
routines in allocating exactly the needed storage (for the polynomial 
coefficients for example). 

6) Isolation of machine dependent operation; most of the operations 
expected to be machine dependent were isolated into small, short functions. 
These are mainly the bit-picking operations and are to be tailored by each 
installation to their particular hardware. By using the pragna INLINE, these 
operations can be made run time efficient. However, since the INLINE pragma 
is not implemented in the Ada/ED version 11.4, it was not possible to 
demonstrate explicitly here that this works well. We were able to demonstrate 
that this is not a difficult or unnatural way to proceed. 

As mentioned earlier we were able to demonstrate that the elementary 
function package is syntactically correct and have been able to execute a 
reduced package containing just the EXP function. Figures V.B-I and II show 


error curves for the EXP function. (These error curves were generated by 


3x with (e*)3 for 50 values of x. Note 


3 


comparing e 
that A= e?*/(e*) - 1 can be up to four times the error tolerance since 
cubing e* sill triple the relative error and the error in @ 3x ean be in the 


opposite direction.) It would be desirable to run several more such curves at 


different (higher) accuracy requests before certifying the correctness of the EXP 
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Figure V.8-I 
An error curve for the EXP function. 


%. REAL'DIGITS @ 5 here giving an error 


The measure of the relative error, 4 © EXP(3X)/EXP(X)**3 - 1, 


is plotted here 
tolerance of 10 
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function; and of course sigilar tests gust be oun on the reat of tne }:orary 
functions vefore they can Se ~3eq with any confidence. 
V.C. Unresolved Prodiens 

There is one major wnacived proslem aasttci ated with Une particu. ar 
method of constructing function pacaages a5 giver. nere. Thia mas te Gt with 
memory space ang tne fact tnat Ate Goes 152 BeeR fo Mave anyiMing a:Brcar te 
tne OVERLAY feature of FORTRAN. [Tn setting sD the elementary Functions 318, 
COS, EXP, etc., & Set of TNedyptes Comatants ave intredute¢e far each Fancti an, 
these are ased tO telcalate “ise coefficients <f tne Pawel pefliéa sf $4 yteamial 
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system has some kind of automatic paging arrangeten:. 
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i ADAfLile: Fumi aca 
5 AISf‘ .e3 Fun als 
LISfiises Furi. uls 
,| 
j oe Yiscletreree "TOUS (aCe ate cl @urmetiors Sh of uevesiaates 
2 od @arly versior oF 3 As Attereroes toy extcerirertas se raiv , . * 
3 ee 
4 badlied Vary cf "me coeste sClS iGe* 1h 8d Face ale wt Oat ewnvard at * 
S oo S°VWC Fal 2 sPOfe Keceseer. €* wwerecose te Fieve Jf fre 
6 oo AQasHl Cn ss iter seraetor ti.t, tre orty ore avalrarie at tre 
; oe Tire © .@0$@ seco Xieg oPfe@ reirsy .fesered, 
' f ee 
t 9 on la -:¢ d 
' tu oe S ESE HO ONE we Pt taet tag oh I%77G7~45 Ft 
11 oo iol a Ce eee 2h ee See oe 
1200 ee 
i3 oo ~ Dt LEO ad tes 
14 ae 1. emer yc at AOL La, Cot m eC OTE “Ay CL '@.e Fre Bare tA Le 
ie) ee fe i am” * ca ee Phe er ee ee 
| 16 oe r, rr ee. ad he ee ee hee eee oe er ok ee Oe et Cal 
t 17 ne 
| lk oo 4.5 Vbeos 
! 19 =e i SORE ON Poe Ea aE Big ee Le geir-t - .* of *¢ valve a 
20 oe le Boe eh FEL MCAT MH TEE ar CH FQ AEN SS Cermart greg 73 
21 oe d. “2-H T Sw, CTF ware 2h AT OTST OrE CF CUM etyemeg 8 
22 o- Be es os THe Fa FO a 8. ti°e@ erprrr *ait, a ok. ar 
23 oe ES I ERS BOE OB PEO roe 
24 oo 
2° 
20 SOO TiC tse. we wet Qs @s'%e. ep? 
27? fcre wt? ce- SE PES. ee 
78 
yam) 
yo bOaCKAIM - Le oP 0% 
>} eg cee eS '+, -°fac(i *s + par 4@ €9Y% nF - + wi 3 
32 Ciree gr: Ferrers fg je ews c+ 8. ) TeOturr BtBl ys 
33 Pinere gr Thee sg ESdt oe thy peepee wean 
zk) PSOE fee POO ee: Fe eee Sep -e28 retor~ ©. ars 
3% einer Pee nae avs bsbo te ey rete ¢{ aT? 
Jo Be My eg LES ee Oe « te aah a®** “at TCL Th FT, eis 
37 ae 1 OD i ee all Oe ee Pe ae ee ee dk refs rs Fo Lb ae 
3” f or cet wv ro se egw lf (eters | retcrn Titers t ? 
sv EWAER IS Tee ee PS Ae wpe refur fF T+ Ceo 
an EIRENE 0s ee oe ee RS, ye PSE retors terry 
314 CURT Pt Pe et ee we A Eee e eet pe vi Qe 
42 PACU or Se See Lee ep ry pete Sere ek) retour: ttt teey 
43 ir hec-d a ire crs or ce! ee ee 2 reterr ee bas 
44 Cnet gor 9 Usb .se her) Pere g t chaos 
a5 FINCL Ins (ts ee Ay fetus ~s \t ge 
4h guner par Bre- "ER, igep sy) retcrr «tals 
47 Curetgiorn Fb. 6 4 ego ae} cetctn weal: 
ad Pryce yer 8 Fe age- alg Frage tat) retitn seai? 
ay Neh bow ch gt th eu egeBayg es Reiger apy retirr eerie 
59 ene tio? wy we EE RSH OE Sea ee ALY retcrn we hye 
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52 
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, YU Adesevl 11.4(7/5063) Inu 26 OCT BL BUstes2e PALF 2 
$3 Sunction cePC Esher) Fetcrr LAUT 
$4 CsMACCAIGAN LatrreFar) return HeAL3 
: $$ CIMELCION Sue li yiwe bt) return REAL 

36 CYRACCL OM ATAL( BR 2MEAL) Feturn SEALs 

$7 Curction SICasheaL) tfeture, meal? 

Se function COS(tsnear) return mEAL) 

+9 function Tas(ssweAt } return MEAL) 

e Cunction aSi®(yr@eat.) return #kaLls 
©} funmecion aCStuseraud return “eaLs 

@2 PuMPction afhrtRgnEal) return eats 
oe) function shvr(eeRCar) return bFAL?) 

64 functtenr ClSme(asnEat) return Peas 

05 CunetCion PFaAvelagagal? return ert acy; 
ee Sanction AST WCUSHEAL) retcen wFaly 
e? CMEC ior ACE Srv sQeAL) return be bl? 
| on 

oe” end eb ruse 

74 

a 

r2 

a | PALO RIOE MF Fae css ay 

Ve 

ry 


Pe ee © 18 overicatet tor tee ene 
>? eo 4 $e sverlioases for Bt 
de or enere £ US ESP Gth and Ww EB HEEL OF OLUAT of BITC Cet ar, 


La) 

ao 

a} F section eC Las iferbmese esavo at) reture “1 ab Ig eof 
44 eo efaits Ester? 

a) ol vs Se relate pe LeUF vers 

ae ent *a*; 

ws 

"6 

a? function 7% Cathe sis Cael chw) ceternr e+ bf Jg oo 
44 eo pra @ JubLTseys 

e? vest FeC ath 4f/we ai. C trys 

Le ex? *s*? 

| 

0? 

9) funceton FOP inl Geese ~sFLUAT) return LUAT Ig eof 
as co rraetae FSP we. 

os eeain feturn FLOAT 1 a3 

4 ens 7a°%3 

ie 

38 

39 function S7*(CasELE IME? ESP stub kh) retnrrs PLUAT 18 oo 
bo *- ofatesa LVULVES 


| EU | nesin returr \/FEVAT CE DS 
102 ena’ e7%3 


{73 

104 

195 finetion “@* (eho Pestle sy «tru cFtcayy return LONSce lat 16 ee 
$06 e- orasma Pel tree 


10} OP aN return EV GchtOarcrypenmys 
108 enti *@*: 
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function */* CALLING FL DAS; TtisTFGEk) return LONGIFLOAT {ig eo 
== Oraygma INLD 3 

0@ 314 return A/SUOAG FI OATC LI): 

ens */*3 


CURCLION SASho2oe AMI CASEEAL) return I: 1EGER ig ee 
"7 tras ge Yelle: 
er a funetion trust suives tne case 7 exit orert 
“7 Of A, SUCNM BTNAL 7,5 <s ABSOX) S79 Sr ASE LZLEXPIN( KL) <€ 1,6, 
ee (orovi4des A 38 feet «D3 
er 10 SPOUT He sacl it? 1ependent ane sritten in essernly code, 
“= pracrta Leld tes 
Zsweol ge ar (Uys 
Jel'.feVwh&~ gz 3 


oesair if 2 /B et eran 
“7310 7 >8 Toy 
1S or >’ ge Z/a? 


ee eR aE 


s 
e dle < 1,8 
len;. 7 $8 o%¢? 
Jog2 © 23 
@nr 4 Leys 
era ifs 
retoren Js 


OND er Ase char tt: 3 


Tarscorstant [ f- ek 33 25 22 PONe: FH ECAYUSe OF 409 in .¥'t coenzler 


functiog. | oe Pee i an eae W's cre) veqcire Lo fh Oss <> me 
ef OF Ain4 bab, dor 3 

7PM AD P4UNCCINY OT ssviler td CH Le)’ eens 

2m At eOakt re otc st 39 asseriv corm, 

ne vir reture F- '8O ts 

Ont Te lPMaernwdop ef: 


function Dede hire me cme aC rst te sch) petope meet {$3 ee 
2 F437? Tru] -r3 

ww A LUPCCIOF ME LeAlLert to nt ratie: OO, Qe 

on tf eo) } we yoenesgr fii aS$e@°r le Cet. 

Heyir return te Al Cr “ee ryp? 


OND Terai. sh ore oS 3 


Function CHepiSe cee ast eesee ald patup, tat {e on 

O2 OT 44 Lsus ies 

P2 A FUNCTION @culvalernt te RITA Ci ree ee AL CEASe oe ak OGM) 
an it edi 4A ne MIpPMmR ect in Q@S$se-iav cr i. a 

oOe@ygin ret itn Rfamm tt EE. FCOUNAGe Coot as. (73)? 

ONS Chew Se oee as 


versit 


FOAM ODOM ORME SORES EONS OOOOH ESHOOE OH OOO BOOS MEOH Case cena eraenavoce 
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lol oo 
Lok ee UfTtiirgess 
Ln ee Sore Of These could «est be writren is assembly code, 


17 == 
tit aaa naa eed dekh ttt ttle ttt tt ttt Te 
L72 
173 
174 funcr {fon Ci rVsig 
t74 HEC LEA Ve 
176 A,YShERL) 
$7) rerorn vArat is se 
17t “e pragna Pei. Iok: 
17% ve gin it * tnen return Ks 
{Ry else reaturt Ye 
| and fs 
ld2 end Cis 
Ie 3 
144 
Bay tunetion cr..f ¢ 
136 Reaithe avs 
{67 A, tS] efhi,eR) 
18s return ["tkGkEX fs o- 


129 \= pranv, Iviblves 
190 regqin 1f aA then return XxX; 


191 else return Y3 
192 end ifs 

193 @nd Chiwie3 

144 

195 


196 function YONCLEANTEGER) ceturn HOOLEAN ig ee 
197 \- pratwha Lyle; 

198 ee Fetter in assertoly cade, 

199 beyin return - mon 2 /z2 Q3 

200 Qn nus 


203 function ROUNNQASKEAL) return INTEGER ig ee 
204 22 pragma Lilie; 
205 ec NUTES Version 11.4 of complLier TRUIECATES TOWAKOS ZERO 


206 on instead of rounding for tUnTrGek, 

207 beivin if KX > 0,0 tren return INTECER(A + 0,5)3 
208 else return INTEGER(X © 0,5)3 

299 end {£3 

210 end’ RIND: 

211 

212 


213 function “*rem"(A, RSRFAL) return RFAL ig e« 

214 “= oraqma LALINE? 

215 es Fetter in assemply {¢ tnere is a machine 

216 -- code tnstruction for remainder divide of FLUATS, 
217 begin Teturn(A © BeRFALCROUND(A/B)))?3 

2ia ena "rem"; 


221 function SIGN(Z, XtkFAL) return REAI. {ig ee 
222 “= Orayma INLINE? 
223 oeain if Z2 < 0,0 then return -X3 
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223 elsif Z2> 0,9 tren return x3 
225 else return 0,v3 
226 ena jf? 
227 end SIGN; 
223 
229 
230 function POLYCXSREA]. 3; AtACREAL) return REAL 1§ ee« 
231 #* pragua INLINE? == 7? 
232 TRAPSREAL s= ACA*°LAST)3 
233 oe jin for J in reverse 0,,A°LASTet 
234 loop TEMP gs ACJ) ¢ X#TEMP?3 
235 end loop? 
236 return TEMP: 
237 end POLYS3 
2398 
239 
240 function JUDIPOLYCXSREFAL? ASAUREAL) return REAL 1s e2 
241 “2 pragma INLINE? 
242 begin return X*®POLY(X*®X, A)? 
243 end ANNIPOLY?: 
244 
245 
246 function eVENCPOCLIYCXSREAL? ASAWKEAL) return REAL is «= 
247 \e2 prayma INLINE? 
218 begin return POLy(X*x, Ads 
249 end EVEYVJPULY? 
259 
251 
252 
253 OS SO OER EE 2S 4S ORES FO HSTESSTHRASS SSHSSTEPSERSSSHSHHReVTVSSERFSOSOUVURBannaan 
254 aw 
255 oe Tne Elerentary Functions itn ADA 
256 on 
257 oe In the past tthe elementary functions have been written in 
258 ©e- assemonly code and incorporated into the niqner order Language 
259 “= orolrar oy peeing Intrinsics of the Lanquage and were nandled as 
260 “= pert of tne corpiler., The acchracy was fixed at the accuracy of 
261 e~- tne nacnhine for «hich the compilation was done. in Ada, however, 
262 “= neitner the racnine nor the accuracy is known until compile time, 
263 ballad [n most cases these functions are napped anto some reduced range 
264 e- and then aprroxitated over tnat range py a Cnebysnev polynomial 
265 we -+terived power series, (The advantadre of the Chebyshev polynomials 
2b6 se are that trey yive a smaller maximur error over the approximation 
267 "“* range than any other polynomial of the same jeqree,) 
266 -~ Tne Chenvsnev colynomials thenselves are not used but the sums 
269 ee of tnetr coefficients of like powers of the variable are used ag 
27 \e- tte constants tn tne power series, Since the Cnebyshev polynomials 
271i e=- contain all %Oners ur to the order of the oolynomial, adding one 
272 e=2 wore tern (tc vet sore accuracy for example) changes all the terms 
273 e- {nn tne nower seres, This teans that the number and values of the 
274 ee coetticlents ust re caiculated at compile time (or later) when 
275 ee reyuutred accuracy is Known, 
276 add 
277 oo 
279 
2340 
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281 
282 
283 
284 
285 
286 
267 
288 
289 
290 
291 
292 
293 
294 
295 
296 
297 
298 
299 
300 
305 
302 
303 
304 
305 
306 
307 
308 
3uy¥ 
310 
311 
312 
313 
314 
315 
31b 
317 
318 
319 
320 
321 


subtype INT.1i.2 i156 INTEGER range 1..23 

subtype INT.O11 is INTEGER range 0..1; 

KMAXsconstant $32 103 

type ALONGIFLOAT {is arrayCIwTEGER range <>) of UINGLELGAT; 


PI zconstant s= 3.14159.265835.89793.236460.264333 
SORT.2 sconstant s3 1.,41421,35023 7309S C4880 160873 
UNW2 sconstant $= 0,09314.71805.59945.3094) 723213 
oe CHEBYSrnV FXPANSIOUN 

2e- Let 

oo £(x) = Sum(4(K)*TCk,x/n)), @-A<2x<sd 

oe C(O¢sk<an) 

ee where 

== T{(«,x/a) = Sum(C Cee 4)#Cx/0) #4) 

ee (U<sf<csKk) 

oa 90K) = C2/pI)F#Intecral(r(xet (kw ,x/ay/Sartcle€Cx/d)ee2)), k>d 
eo (el1<ex/di¢zl ) 


C2/rtv*inteagrall(ttateesly) )*cos(« Fy) ) 


ee (O<saydcer1) 

oe t(x) = Surf Sum(e (we JCC ,pv*CKSA) BEF) ) 
oo (o<sk<¢sr) (U€szyerK, 

ae 3 SMCACV) x44) 

eam (Uc=sjczn) 

-_> & anere 

oo ACH) S (14a) ##44+SyP(R(KIECOd, Id) 

ae ( )<sk sn) 


Packaue FFALCIU fs new FLIATOLC UC Ke Al Ds 
use RFALLCLUS 


procedure iCAcSti4¢ 
Asout AMF AL? 
MAb fw cebe'Qrs 
Fe a a eee ee 
Ce ENT Oe). 3 


oe ly 


1 of 
Po 
« 


% A3ves even ;ower series, 
} Qives oO}. power series, 
m gtves tor even and o3d4 terms, 


oa hy 


s 
I 
i 
{ 


oe t, 


a 
— 
= 
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338 
339 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 


\e= Le#tl, I 21 is undefined, 

ew Tne CHEHYSHEV constants: 

function CCN, JSINTEGER) return INTEGER ts <= 
“= This fs not efficient, time wise, put no matter? 
“= only used to initialize package constants, 
eo never used in real Cime, 
e- Tnis IS efficient space wise, since it eliminates 
== the need of storing the sparse matrix C(N, J), 

peyin if \§\ < 0 or else J < 0 or else J > N 
Or else JDD(neJ) then return 03 
elsif (Jsu and N=0) or else (Nai and J=i) tnen return i; 
else return 28C(N91,Jel) © C(N@2, J33 
end if; 

eni C3 

TEMPS LUNGIFLOAT? 

CTEMPS INTEGER? 

begin ee TCH.S!I™ 
for J in ©,,A°LAST 


looo TEMP s= 0,03 
for * in J,.,A°LASS 
loop CfreMP g=2 CCL#&he], EFTeT) F 


e- Version 11.4 wont multyplyv py floatina zero, 
{€ CTkHP /= uv then 
TE“P $= TEMP ¢ CTEFPFbE(A)? 
end if; 
end loop; 
ACJ) $= REAULCTEME ER (1D ee (CML eJeo] )}3 
PUTCAC I), alOtTH => 15)3 
end loop; 
newline; 
end TCHiSJU3 


funertion LI*(BsALLOXGOFLIATS GsSREAL £2 10,0) return IvwTkitk ig «= 
KSINTREGER 32 13 
KREALSDIGITSSINTESER sz INTEGCEKCREALTDIGIIS)? 
neyin 
PUT(Ca, WIOTH => 15)3 newrline? 
PUTC2Z4#u, wIUTHn => {5} 37 newoline?s 
looo {ft ArSCREALCBCK))) < V¥EFCPREALCNIGITS) then 
return aA © 13 
elsif * 2 R°LAST then return Ks 
end if? 
K $2 F @ iy 
eni Loop; 
end LI“3 


ewe Hote Boxxx(e!) = nagnitude ot range of varianle in fiunetion xxx; 


as ie 

ow HesIN( =f} = pis4 RACUS Cet) zbei/se 
oe BeASIL (1) z= 0,375 FoTAUN COL) z o1/u 
oo beSsuRT( od ) = 0,25 MHeEXP( 1) = 0,25 
ae DoATAN (Lt ) = $3rt(?) @- 1 

= delbN(=) ) 2 (2 = sare (7))/07 © surt(2)) 

oe Sn5INn( 1) = 1.0 Ae cUSK( 2] ) = 1.9 

oe HeTAVHi( ot) z.9,.5 roATANH(@} ) = ).25 
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395 o- teASINH( 1) = 0,375 
396 
397 - 
398 2w2e Tnese constants are telieved to te eccurate to the 24th decisal digit, 
399 
400 GeSIMNs constant AQLUNGLFLIAT COL, KMAKR) 3@ ( we Initializes: 
401 PI/4, 0,.94977.049415.68744.77636.82691, 
402 “J ,049H 4.0412 9037030065401049298, 0,00038.077334.36182.82730.90288, 
403 “0 .,09000514 395089091293 Zug. 90335, 0.94¥00009030.736812185446.85472, 
404 “0 ,0000U c0UDU0LO8 I 1b 3059782291, 0.9V40000000,00004.27564.99507, 
405 “0. AyvVNHU. J0VVULONUIN L903 436072, §,.00000.0900000060,00000.17640, 
400 “UU N00UDLVGUIILNOONGO LL CODULONNNS, O03 
407 
408 BCOSsconstant AabLuUrSceLIAT Cel eahax) 32 ( ef Initlalizes 
409 PI/4d, Ue9dJO32I9LITLIOF608 00127004060, 
410 MOL 4ASaLOOAEILIO FH LhOII2LUTIIAING, 0500139 20144930113814064079.69979, 
411 @0 ,DVU00N YIN TYLON HIN 29390 L0UA/, 4675304 000275276595260 11873952, 
412 MU -IIUUU LOI IAU LC a7 IDI IGKOBLLO4dada, %.Y¥0000200000000055.4699548.54157, 
413 0 VOVVIN LP VIDICIWMAGLISETIGLIIVIVIE, VU SYVLENONLVONDILOUNANULIONIOILG2990, 
414 PO, ODI LOIN IO LOU ING Vata OSs, Oe): 
415 
416 Be FANE CONSCANE Aged SFL IAVCeL, eK AR) SS (€ om Initialszes 
417 PI/t, LeH2ZHYSIHAOSSCIOAIIATSI3IL GL 77491, 
41+ OU 2TIALMUSHIL INO Ma CAREY SZ UIQIH, Pe IOV LG FHA 4447 LH 90730239490, 
419 DeQUNAIO LMI LOST ISATLVIAT I Coad dy, CIV ONQLILE 2429567027 961097971, 
420 VeMVVIG STIG PAD IIOL SSI HAIRS S, UM UGADLUOUMILAGSAZLI4IIV2AT IST, 
421 Oe O0UELCHMUIO Ves TIS ZIT Le, CEs GONCAVIDU LUC NODC oMs3zTO709 94, 
422 Oued crepes me CIAIIT CAA AD, eC HME DO COOAOD GL CHUA UU 200 Cte 2a) 
$23 
424 BeoATAsconstanl Abe Gor UCAT Cele er vei ck) 8 ( of Inittalizes 
425 Vet 42 ad oO23 C7 uM dr NL ERAT, re PAST CADIZ LAW AI LL ANA SHLAh7, 
426 U0 ,0 25999 L 27204 Cn dh deh CVI T5, He POECRAS 29 THY QZ LOY I2 HOGG 32, 
42] SUS TIPIP SII GCM IIIS LONI Liao, TUIMDNLOSAIS CT nraS¢orUloteadsggant, 
42 SQ MIDPIVDE VID PSL Kah ISA LSM IDS, SetE STII LOOVIS KT IZ REAM DL 2ITS1, 
429 NNnsuALOH Pu bes IL ATL I Id SG, Co IVC IGG ACO I PSS AQQ2oIUC8G, 
43vu WO SVQ PaaS to 2757, TIC IFN ILININO DV Qe aM s2TC IRSH97 ) 
$31 
432 BOASTASCUNSCAaANE Veu e er tAT (Cele et 82) Ge ( we Tritializes 
433 1,375, Pei Us bao dveehhstanlreowg3e i, 
434 ToVEZ s0o 2A! 4a rd PAT HLS Jr, ei MID POYS2H toners ele dIILRIIPGD, 
435 VeliVOUIL ET GOmMle steve lp slsavel, re ETE AT I Ce sal ac sSohoa lL ISi3, 
43 VeOI sr suurr Zs Feed Css sulce Hass, tag dD DOM VO STS & Jah gael i, 
437 QAI IOLO MNO Le PIP gs hoe seu, te HCO LAID Tee me PES ae ?, 
438 OPPO HVLINO ILE PIL I ealeed at, pe tO VOTIVE Fe de PV AR ZL IG271) 
43% 
44.) HosbwesCONMslant woh! ahi SAT LOD oh tt) 8S € om Foultfalizes 
44) 1. Per MP SQLP ITAL ZS ISS LaF: ILIGSEG, 
447 OLDAISGDLASAVYS EAD ICA7ITIS AOU, Oe eeu ILIA TL se Moise IS eeT, 
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444 


n45 function fArn(xrseFalt,) return REAL Ig o@ 
844 22 prayra Itulyr3: 22 ? 


847 Zieh aAL $2 xX 

R45 oeuln loO>p 

84) tf AKS(7) <= S*PI1/8 tten returh TAPLSMALLCLZ)?3 

aso elsif &4S(Z) <= #1 tner 2 t@ SIGN(Z, ABS(Z) = PI)3 
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ee eae ee on ee e/ CT af Lornes Ze rAGFE 4 
©. .- fec 1 78 Jy? 
“~ a 
ed a - 
oe . Yuaf ' Base 2 tan -§ 
oe ™syt ¢ ‘ * («ec @ SS jP 
oo neat bet S Sait @ Foe rare @ 2) 
oe airy? s t @ ee a oe BEE Ce rgd 
22 St ise” e,* ', rr af eae, rj fs 
en - 3 } é ] “ % 6 ] “ 
oo mw, thiverr) = 4 0 Sy¢ Pe! »,! ''.5 32.3 24.2 is. 
ew on On: oa’ 
@ee@ eveaeusndsewesweeweeeaews et @®@ae@@eegeed¢dgeaew#s¢e @eeg * @ae w2@aeneeeeeaeseeeeaedeueseesFsregeeeseeewvweua@eeeeeeec&ke @ @ 
t yo Atae th eerve at ) vreture erthj. Je of 
e uratre | { -3; #2 ? 
S77 At gS avi ¥)3 
cecvt. t@ ? 2 er Fe sce vasy Chey 
retrr tobe h ea, batlan)?3 
e'stt 7 €2 Ve fi CATA. CL AA® Cher 
Teter rer PUN OR HL ZGe Ot CHOCLYC(Ze1L 00740701 ,0), PLATAS) D3 
ey: "£3 
PMC TN bee FESQ2 @ eM hey) et 47, PLATAR)) 3 
ent naA.? 
oe Cramoutealictay tires 
oe nest? P i * (a > “jn 
oe Ave? Tet A @ of 6 8 e@ CA © HON @ EY?) 
Sd eorset: 7 @ & & (A @ MYCN & 2) 
oo iart, ¢a3@ 127% (ranges, tar ci/se))es 
pe ns 2 3 4 % 5 7 # 9 10 i | 
e\e= eforlicerr) «= 4.7 H.3 Fo? WSF LH.27 32,1 #+$$3.0 #$15.0 16,5 17,9 
eo 8 fw pathos To, Sufttsg) | Lh 
Seer-eneaqgewnwzweewenwe@eaeewnee ease onse@ es Goe@we Bs eosGeaeaeswes es Geese eeseese ana @ese eee ae ee es G@uweeeu Geese @& we w& 
ae Nie dus 
@e#@q@weuweweeeewreese egqoeweet @G@eaenwraeaweeseaeeoeaeaeesaoeewGaoaauw sees eaueeaeaescgesgee ee@SBaesesesae @aecue S@Oeeoaee &@ & @ SS 
fucction ASIA(rsRKEAL) return REAL {8 © 
ee oragrta IVLINES e2 ? 


be 31 ee aseert AtS(1) C= 1,03 
lf ABSCY) < FLASTeLRAN then return ODOAPOLICY, PAASIN) #§ 
eisif Ax5(V¥) P85 2.0 then return SIGNCY, PI/2)3 
er-4 if: 
return ATARCT/SSORTCC1,0 © YROE(1.0 = Y¥)))3 
era AST? 
ow Carrutation tttes 
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NYU Ada/ED 11,907 7/0/61) ras gs MCV ab Peeross,s Biter V/ 
908 co Resets €* @ C4 & YN 
909 oo Aves Pele 25. € USE gn © Teeter © 
910 oo heISES © 42CP Lalas «¢ t feqrt) « 173 + 4) 
9it oe eorst: S09 @ © © CA @ SJE fatasut @¢ anlsercj) « 9) 
912 
913 
914 See ge@e eecdoew eevee eee e Gee ea GeeacneaoseeGoeoeceFowaewswss@ FPeSGe SOOT Ge eGeTUa Ge ee SFaoouGweesee® eeevesv ae 
915 ae 
916 oo ACus 
917 oo 
938 @@eeeewsedseces Bees a2eeaeeesteeeaeoaeaegeveeaeeeuese@eeGova vee FCecees@eewt aeeadeseoea Gaeweeeeese*Zoaeeaeeseea ea & 
919 
920 ACGS Ctr «ST2ReaAbs 
921 
922 tunctton ACUSCi sar Al) reacurn «hal {5 oo 
923 e~ prages FIle; 
924 YtkeAlL $= 1, | U3 
925 begin ee assert Antti) <2 1,3 
926 Lf ALDSCY) < ACUSSCUNST Prer rperurn SORTCY9CZ.9 © ¥43))3 
927 ent if; 
929 return PI/2 © AsteCu)s 
929 end ACGS? 
930 oo Co*mgitation titres 
931 7° ‘ests 2 @ 2A @ (A @ F)(thlsuart) » 3) 
932 == Aves 2.9675 ¢ 2.80258 + 1,298 + 
933 oo VeTS(A & WI(CnNCatean) «¢ nlsaqrt) «¢ nflasinjy/§ © 4) 
934 oo worst: 30 @ JA @ (A © WC nlarar) © anCsuart) ¢« $) 
935 
936 
937 Sse @e@eeeseSBeeseee ee SG eOeeae Gee ee Gea e@eeaeseoaseeee estes G2aeees eee ease Goesereseae ee ee oantoartaan a2 
938 oo 
939 oo Slt 
940 oo 
941 SO SSO CCS OFFEND SOHSFCSKHFSHSS SHSHASSSSSSSSTSSSSTSSPSSTSSSSFSVASSSSSHVHBSEBETGABZA®SZE®S 
942 
943 function SINWCx2RE AL) return REALE {s e# 
944 ** prayma [NLIsty -- ? 
945 EsREAL} 
946 begin $f AnoSOX) <2 BIST WHIRAN then return ODOCPULY(CX, FaSI hr) 3 
94) eni its 
948 E s= EXPCA)? 
949 1f £2 0.9 tren return REAL CREALYS LARGE) /23 
950 end if; 
951 return (& = 1,%/8)/27 
952 end SINH3 
953 
954 =e Hart, page 104 (rangeld,9,5))3 
955 oe yn 2 QO j 2 3 4 § 6 7 
956 ewe ejoylucerr) = 1.7 4.2 7.9 19,1 13,3 16,7 20,3 23.9 
957 oo Rangqe(O, 1,0]3 
959 “oe 3 3.0 5,2 7.7 10,3 13. to. 19,3 
959 
960 
961 Sec e828 0 G@ESeOSGSESeOGSG@seSV eS ees OOS es eear Ses GBS eaGQ Oe ese avsiese@iasaBoeeasQesacdse @G 
962 ae 
963 oo COSH 
964 oe 


a eT ae + ‘ re eer rs 
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965 
966 
967 
96a 
969 
979 
971 
972 
973 
974 
975 
976 
977 
978 
979 
985) 
981 
982 
983 
994 
Jus 
9Bo 
987 
933 
989 
490 
991 
992 
993 
994 
995 
996 
997 
994 
999 
1909 
1001 
1002 
bou3 


SE 2O@eaeVaeen ote FT Tee SOA Geese ass eEAsB@FAsea BSF SE 8 E RSS FE SSSSESHVVSeETSeSESGaSEaoTUse eee swvdoe®@ 


function CUOSH(ASREAL) return KEAL is «<2 
=< prauta FNUANEs ©= 2? 
ESREAL; 
begin L€ AnS(x) <= BOCUSHLRA then return EVENJPOLY(X, PICOSH)3 
end {£3 
BE se EXP(X)3 
{if & = 09,0 tnen return KFALCREALC LARGE) /23 


end if; 
return (£ ¢ 1.0/8)/23 
ena CUS? 
= zange(9, 1.0): 
oe nN i 2 3 4 5 6 7 


oe ejo 91N(err) = 2.3 4.3 6,7 9.3 12.90 14.8 17.8 


oo Talli 


Se se eae Se SSO Os SOesB OG SST SESTSSEPTESEBSSSSESSE EWA SEFSEBTSEST*SSSVSRAVEIVSRBSISSSESSEVSISESSBOSEGsS SG @ 


TAWHICUNST3 REALS 


function TANUCXSkFAL) return REAL fis e= 
ee fraans INLINE S 2 ? 
Bn2XgRrAu? 
yvejir if AnS(s) <= RL TANHCRAN ther return OOM.POLYCX, P.iTANH)? 
CLsif AbS(X) >= TANHCCOVST fhnen return SIGK(X, 1.90)3 
ena if; 
B2X% 32 eral 24x); 
return (FAK © 1,Q)/CKE2X ¢ 1,0)3 
ent? TAr43 


SH eG eeSe eG SSeS FW VW GBFeseswe ses ee woae* Gewese SVS FG OSeSeS esses SEGos een snus 
oe ATANH 
oe 


See GQeweweaeet ase Ga OS Ge OGee@enOOeG2SGor one ss FSG a GQSeSeSResesesQOaesesee 2a et seennesu & 


fuactian APARNCXs&*F AL) return REAL JS oe 

e- Jrasi4 Leb. gs e= ? 

ve jir “- assert AkX®S(¥) < 1,03 
if ARS(K) <2 *lATANMI RAY then return GODCPOLY(CX, PoATANH) 3 
Pisift AMBCA) 92 2,9 then return SIGN(OX, LN GREALLILARGE/S2)3 
evi it: 
return GeCCte& ¢ AYSOLLG © X))/2)3 

Oni APar-4+3 


oo Kandel), 251% 
a | 2 3 4 5 6 7 ] 
= elogyiocerr) = 3,4) Sel Tey 9,7? 11,6 13,4 13.3 17.1 


ee a 
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1022 
1023 
1024 
1025 
1026 
1027 
1028 
1029 
19306 
1931 
1032 
1033 
1034 
16035 
1036 
1037 
1038 
1039 
1040 
1041 
1042 
1043 
1044 
1045 
1046 
1047 
1048 
1049 
1050 
1051 
1052 
1053 
1054 
1055 
1056 
1057 
1058 
1059 
1060 
1061 
1062 
1063 
1064 
1065 
1066 
1067 
1068 
1069 
1070 
1073 
1072 
1073 
1074 
1075 
1076 
1077 
1078 


oe ASINH 


Ses Oe es es eGeasee SSB UeEsee@Bs OE e2eeesV ess Ete BF FH OF GFT FSEETSFTeCSB A Qeoeeoeees eg Fsenn UE GE Ge & 


function ASINH(XSREAL) return KEAL is == 
“* pragma INLINE; == ? 
pegin L£ ABSCX) <= ROCASINHLRAN then return ODN PULY(A, PAASINH)? 
@lsif AKSCA) >= SQRTIREALLLARGE/2 then 
return SIGN(X, LNCARPS(X)) ¢ LNO2)3 


end if; 
return SIGN(X, LNCABS(X) + SGRT(1.0 + XBX)))3 
end ASINH) 
« @ Rangelyu, 0,375)3 
se ns 2 3 4 5 6 7 4 y 19 


em» elogiv(cerr) = 5,5 7.1 8.8 10,3 4.1 $1365 $15.0 160.6 18,1 


= ACNSH 


ACOSHYCUNSTS REAL? 


function ACISACUSKEAL) return KEAL is ee 

ee pragma INLINE: e= 7? 
YSREAL s= YU = 1,03 

vein ee assert >= £.03 
{€ UV <2 1,0 tnen return 0,03 
elsif Y¥Y <= ACUSHALCONST then return SIRTCLFC2.0 = ¥/3))3 
elsif U >= SGURTIKEALCLARGE/2 then return UVCAKS(U)) ¢ ONL2; 
end if; 
return LACSGRTCV#(1.0 + t1)) € UD? 

end AcUsSn} 


THOCPLFIVE3 constant PEAL t= 2.57 |= seeded by Conpjler version 11.4, 


begin oe Initialization of package tutFun, 


TCHSUMCPLSIN, 3eSIN, 2, 1)F 
TCHASHM(PLCOS, teCOS, 7, 0)} 
TCHUCSUMCPLTAS, AITAN, 2, 0)3 
TCHILSUMCPCASIN, OASIS, 2, 1)3 
TCHULSUY(PLATAN, HBoATAN, 2, 23 
TCHISUMCPLEXP, REAP, 1, 9)3 
TCHCSUM(P.SIAH, ALSINH, 2, ODF 
TCHUSUM(PLCOSH, ALCOSI!, 2, 0)3 
TCHUSUM(PLLTANH, SoTANH, 27, O33 
TCHISUM(CPLASIMH, HLASIANGN,2, 0)3 
TCHISUY( PLATA! , BOATANH, 2, O)3 
TCH.SUFM(P.LLN, Reolitt, ee 9); 
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1979 
108u 
1081 
LVd2 
1043 
1044 
1085 
1086 
tus? 
10569 
2RY 
1090 
1991 
1092 
1093 
1994 
1uys 
10%5 
199? 
199R 
1994 
119: 
113 
Tlie 
1103 
1194 
) Bebe) 
Llie 
1197 
Livk 
tiay 
1199 
1121 
JitZ 
$1143 
i oe 
1115 
Lite 
1117 
1ii* 
lalsy 
112") 
1i?71 
L122 
1123 
1124 
13125 
Lido 
1127 
1126 
i 
1130 
{ist 
1132 
1133 
1134 
1135 


PEnHoSUMCPLSUKT, ogSURT, 1, 9)3 


ror iin ©, bit Cs.) 
Loop Poloedcl]) 32 2*rPoL'C))?3 
end loos ; 


ACuS nC SI $= FO PTAULPLF LVEF ee eREAL DIGITS: 
AC ISG CSE $e a. Os TaClLVor EVE eer CALE DIGITS; 
TArtoCusS? 35 RCALC BEALS LICTTS) 981.2 ¢ DNI2? |f 1.2 “= IN(19)/2 


POULTCATHSCCIST, «LOTH => 15)3 
PUTCACHS FLCNEGT, AEPTH => 595)3 
Pep Pas 4oCOUaSt, #«TaTH s> 15)¢3 
neacline?; 


bk P ALLA eh SS DIC KRPAL CRE ALCL ARGE) )3 
SuURPTORPALCLARGE 3s SQRITCREAL CRE ALS LARGE) )3 


Milt Ci CPE AL UAB GS aliyTrl => 1S); 
PUT(SIPPONFALCOAHGE, ALOT => 15)2 
VeV GL IPe7 


end Fle Kits 


Sea ee Geese Geese et ae aa eS aes ee OG @aanGeGsOQG*oeageaaeVl eee aang eeeuseaeae Sees SF Se BESS enwe SEB ee & @ 


-~ ee kerb (rg 

lad nice, Jonn «,, "fue aAnyroxiration of Functions", 

o~ Ridisonwvesley, 3964, GALA) KS VI 

ne dasrairiyus, Ceci) Jr, *Antroxinations tor Luiteal Cornuters”, 
-< Princeton tiniversity Press, 1955, WATA 37 ©,2 

oo Hal stor, Aintrony, "kA KFirsr Course in Numerical Analysis", 
oa sesurasertll, 1IAN, JACST ot 4 

= Sivier, Averv “artinr, "*Cnrervsnev wetnods {nn tunrerical 

-- Approxfoeatior", rPrenticeerall, 1466, QAZ2i .S65 

2° Provanseii, Aleney sikolaeviten, “The Applicatian of Continued 
ee Fractions ant Tnelr venerait{zations to Proslens in 

o- a, praxiration Toeory*, &. “cordPoff, 1963. A221 ,~K473 

oo “art, Jon be, Ot ale, "Computer Anproxinetions," John 

ce 49008, Irving Allen, “tumerf{cal Analysis for Conrputer 

oo Sctence," ‘vortrenollani, 1%7A, Qa297 N59 


Se eee G2 Oe ee asa eG Ose Soe Ges G@aaa Bees eananesgses SBeu_se OeeeEeE ea acewes seen eeu saenuzea ae we 


wit? ELeak os? 


proceisre teiw® Js 


NYU Ada/eb 11.4¢7/96/81) THU 29 OCT RI 10294328 PAGE 71 
1136 type REtL Is divyits 53 
1137 type LONGIFILOOT is digits 7; 
1138 nackage VEwLFLELFUN is new ELELFUNCRFEL, GONGLEECUT) 9 
1139 28@ SEALE LELPUN 
1140 beqin ee TeSf 
1141 PUrC*® TESL OK %)3 
1142 newwline: 
1143 end TEST? 
1144 


No parse errors tetecteu 
Parsing times 732 secon%s 


no semantic errors detecter 
Translation tines 2238 seconas 
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APPENDIX B 


A Listing of the Provrams and Subroutines for 
Computing the Chebyshev Coefficients for the Various Elementary Math Functions 


These programs are in FORTRAN but for documentation purposes, 
the FLECS listings are first piven. The numbers in the right ‘1nd column 


of the FORTRAN listings show the line tc which it corresponds in the FLECS 


listing. 


— se ee 


COFY » OSUBOE 


00001 
09002 
00003 
00004 
00005 
00006 
000? 
00008 
OO009 
OCGLO 
OULI 
CVete 
YOuLS 
wIQ14 
vOO1S 
OVO1s 
OO ITS 
OQYCLS 
WOOLY 
CU0LG 
od 
aye 
p85 
wu, 
OOS 
COIL6 
WOOL s 
VIVES of 
OOOLY 
000 f° 
OHO} 
wuse 
wad 
ree 3g 
yep hs, 
age dé 
Guys. 
yyy de 
Gon ty 
Pt raie 
Ciyete 
ru. 
gas 
yas 
VOOAS 
HMgade 
O04? 
WIPO Tal 
WIL By 
DIST eaeA 
099%1 
gs 
duet 
CHOPS 
HOON. 
CHT yey 


“ee 
hire. ye 


aera ee Pewee ee OK Sew ee ET eee aT ee ee wee eee eee ee oe 


OO OTA aAWIA AUG AO 


woe 


~ 


Lo 


THE FROGRAM TRIGOE CALCULATES THE CHEBYSHEYV AND POLYNOMIAL 
COEFFICIENTS FOR Off OR EVEN CMOSTLY TRIG) FUNCTIONS. 

ALL COMPUTATIONS ARE IN FOURLE FRECISION EXCEPT FOR THE 
ERROR TOLERANCE WAICH IS [a Olwole. 

THE ACTUAL COMPUTATIONS sRE CARRIED GUT ln SUBROUTINE FOE 
THE MAIN FRUGRAM ONLY LORS THE INTERACTIVE BIALOG TC 
DETERMINE WHAT 15 (0 tt CCMFUTED AML! HOW ACUURATECT. 


THE KANGE UF THE FRROY LTO ATION Tl THE FUNCTION £5 FRO 
-at ft AR, ote Tat USER SFEC ir Te > Ab =. Ur LESS THEN 
A DEFAULT Ab wv $5 Ucth. 


g 
FROGFEAM TRIGOE I NFUT SOUTFUT | g 
TMPLECIT Guvbhe eRe ii fd  A-Hee-2) z= ¥ 
COMMUN We abe | -& 
EXTERNAL FoiNts 2 lide : 8 
FT ANP: ccc etd me x 
bye ENB Der Dat Ass 
FAL oR i oa o 
ETANH' s TAs a 
Dang 
FAS| hes atitee | 
FAI + tte ne 
relist» 1 {%% 6 23 
FaflANls rive» UZ 
bALNT» ral ies 
PS bel oe Stdeel 
Fite > a1 NHe 
Feveat oN 
PRE Ghee Rad 
ete PPS  F8h ee Na Seay Bel Cole 


WHIvF 6. FUE 

PRENT @» "ua RMEN -LS PpGe ae 
. Ereafi @e | 
» FREIND 199 J 

FORMAT yds 

FRINT 8» ° 15" *. 

Reef &e at 

RRIV) @e oi 

en? | See eee ST 

ri CW. ) das iNet 

eo oe. 4. SR) alte Boots. wih * Se SHO 


rey Ct ans roe WHT Pe ae | 


eS wt Seok se 
’ bua? 1 we le} 
9 ee eee re ee hehe Bere 
> ° ‘ yer yey eo ras fhe 


’ t) re [ft 
: ; (Yo ,Fie 7 OH dH 


. ° ; Tt seek Fb. Geet oa eo 
‘ © Lebi bot cr TANMS 0 ABS 
veoh THe 
Pot, Pee Ten 
| rr at rT Pee Le 
Ceth FU F cfd ig PUN Ge 
: poe EN 


00083 » « f .€0. 10#COSH 

0OvS? » 6 « YF AR GLE. 2. of wn = 9.900 
00060 » « « CALL FOE FR Cirats the 

Q0061 “< 6° Savin 

JU86. » « "fF 6G, 1GRS0RL 

ON0e. — 6 YR CAR wth. o. tee e 25 Ae 
Q0064 CALL FOECF SUA *s SOA. 

09344 » ». $-FIN 


DQ ee 
e 
@ 


Qyuues » oe FLT CER. POM] Ds 
VIUe ” a? Sr) rd O rs ec | 2a, | 
yVod re © \OC a?) a eee, Bee 
VIVoY ye: ce tee LN 
aye ry ; » lf EG. LOMAS in 
acer Tt g. os a) ee ee eS red Ce 
Cys apes »~ 0 6 Come POR Fa teie Aah 
i | eae a (4 ee N 
ou 8 ey kd, fmt et 
elke SS, ss : ; if a re ee ese 
Pee ~ oe lw Ltt OR UR CR TAN Toate 
Vow” ae ee 
ny 9 » 6 | sbte, J uedee Toone 
nd 2 -o & “Pe OAR she ofa wl Sak TOL. Ghd: - , 
ere re er See} ere ey 
7) Le : vo eb Uh 
eS 5. 1 wba. VO MRes [te 
( rs Sy) ee PP CY Se uhes 
4 ‘ ‘ GS) ee 2 ke, Sey oes 
of yf 
; os eo eh th, ol 
ts ah: Og i‘ eg 3 
: os, RMR Sct Ts im 
, is 
a ’ , Aad 9! Oz 
t sae Ae Sear SG ‘ 
pescab TO ie : oer eee ee 
ney Mee or Bb a te ce te 
ng fh23 : ' oe & Le 
hooey » oe eT ESE OR EMT Be TP Lk mre adie Os 
» ote 


Vuh 


Bre ERS Tite ~ 4 


i 7° ¢ ns ra rt Of, Ce Ce Eo ‘ih OR oe VEN UME FY CHE Y eatali 
; : the yy oe es ES | at SOT STEL Piuec Tob. 


. i$ -\ fae ' parte ° . eer 3 a meg fi - if s]N. (i tye reir | 


oa, Poe OE . ~ 43° ~ $ ). 4 1258: iS 6968 5 HE os yi? 
ms : cb Este tobe tht Ug & Tee te eee GEE ced TRE ME 
ae DM wee HES ce ipbbe fetes. rl’) Ue, rt ee arco) ee. ee TaMiin. 3 
hoop Mp eT eR ERG ENTS BEG Cee RT ls MER tte 
ae a te 2, at ey, she PENT eg ok NE ARE 


iT od Se a 2 Mw fie, 


met \ ee re So i aa {she é iC Tg ats th Me 


T 


ri, aes | oe 


OOL14 
00135 
OOLL6 
00417 
Ov118 
VULLY 
0010 
Qybtl 


AS [oS Tre mbRar OF COEFF EC: 5's FOR TME FOL TNOMIAL 
AbGKOK (MAT [UM Ub 555a-0 COMute., iN ACOEF. 


Tre Of fume lus «SiMe Tame Ss ATAN? ARE LC OWMERTED 
Th eed FONT Pe stm be Ce a he ASP: Ke AML 
tar a) 8 Ie ADEA Pal bat eh at lve ERROK BE UNIFORM 
peste ibook i et “pee Loe bt pee, Ahh ALMA tie, 


Me SM ge Pla amet ee ON Sa wa 


WIR. 
ywohrs Tre ke LE a bt boa AAT OM OS VERIFIED 
vs or oe Boose eb ees fe. cele tem me Om OME fhe eva G Py ae. 
OES & 
DOTA SUR RE Tea ee ew ths 
2 aa PAPEL 2. Celtic oh) “SEN a ee Owe. UE 
eS {PUPS fee teh oo 
PRLS DEPeM fre te gD 
vi bao boaly wei Ms - ie ® es ae 
ea 
Ly id Bee, ee en ee 
trol ff ir nS | Pio an Ce. lew By 
wel da 
Nee Me : 
rey 4s oe Te ey ee r 
wes rn or ‘ 
St fd ee a 
4 
1 . - ® 
ae anes a 
ex Las t 
.ae 1h oe ie ~ i 8, “t. 
re ae r _ > @ORGRSEEReed 
4 SS eck oe “a, - 4 i 2 © ee © Pe 2 Le 
3 
a 
‘ wre pf wot & . 
e ‘ 
» % 8, ¢ 
¢ Pr: fe» * a oS van 
mh G. 
ee 
\ ‘ , 1% ~- : 
Secs. ; 
' ' . ‘ 
( ‘, 
: ate it % 6 
' * q 
, ees! ph ve °55557, 8, FI 
' : ' fy Bg - fhe 
re i + ‘be . 
eee rim Me, ‘ 
i 4 
1 Ny aa 
' oA Ae 
ye ay - ft otee _ gas en © ee 


my ea > [Na @.- Me . *., : ber a La! 


~~ 


~—* 


‘y Bs 
, 
a 
¢ 
“ 
rT] 
a 
“ 
7 
» - &, 
- 
“4 
ay 
Fy 
mw. 
~ . 
why 
an 
ys . 
. ’ 
r ‘ 


’ * a 
‘¢ ° 7 
. * 
e 
1 
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0007¢ 
v007. 
9007¢ 
0007% 
Qddey 
HOE} 
08087 
003 - 
99084 
00085 
00086 
00087 


wOOSE C 


00089 
90090 
0G! 
00092 
00093 
JOG 4 
Y0075 
0009¢ 
My? 
a 
GUO 
OG idY 
oQi. 

2010. 
00103 
00104 


Q010% 
00106 


END 
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C 
C 
C 8 THIS FUNCTION DOES KROMBERG QUADRATURE ON THE FUNCTION F 
C a FROM THE LOWER LIMIT A TO THE UPPER LIMIT & AND TO THE 
C | REQUESTED ACCURACY OE. 
C 
FUNCTION ROMBER (Ar Be Fs OE) 
IMPLICIT DOUBLE FRECISION(A-HeF-Z) 
U GOLMENSTON O66 20) 


INTEGER FOUR(21)+5 0 
GIMENSION Q620)+ Q0(2) 
EQUIVALENCE (00(2)»0(})) 
C 
UATA O/0/+ LMAX/1/s FQUR(1)/4/ 
C 
Q(0) = (F(A) + F(B))/2 
H=B8 -A 
IMAX = 1 
C 
BO (l=1920) 
> IF (tL.GE.LMAX) 
LMAX = L + 1 
© » FOURCLE#L) = 48FQURCL) 
e eeetl IN 
C 
, MEHL 
SUM = 0 
110 (]-19TMAX) SUM = SUM 4+ FCA + Hec2e@] - 1)) 
- SUP = SUM + JMAXRQ:L-1) 
» ]Ahx = 2RIMAX 
Q¢L) = SUM/IMAX 
lu 
DOCM=1 rb > 
N=z=i- # 
©» « Q(K) = CPOURCM)SQC(KA4L) - QCK)ISCFOURCM) - £) 
® ee of IN 
C 
ROMBER = (& - A)KQ(0) 
C » TO (h=ieb> OF(K) = (Q¢(0) - QUK))B(B - A) 
. » FRINT®, GOMBERs (OUCA) + Kelel) 
¢ ‘ 
» O0Q = DABSC CQO) - ACL) CB - AD) 
© IF (ODQ.LE.0E) RETURN 
i, i 
oaet IN 
KE TURN 
ENL 
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00113 
00114 
00115 
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QO117 
00118 
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00120 
00121 
00122 
00123 
00124 
00125 
00126 


00107 
00128 
00129 
00130 
00131 
001352 
00133 


00154 
00135 
00136 
00137 
00138 
00139 
00140 


OO141 
00142 
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FUNCTION SIN(X) 

IMPLICIT DOUBLE PRECISION(A-HrP-2) 
SIN = DSIN(X) 

PRINTS: °° X =%s Ke ° SEWCX) =" SIN 
RE TURN 

END 
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FUNCTION TAN(X) 

IMFLICIT DOUBLE FRECISION(A-Hef -2) 
TAN = [SING Ks /ECOS(X) 

FRINT@Se ° © =%. Xp ° TAN(3) =%s TAN 
KE TURN 

ENG 
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FUNCTIGN ASINGK) 

IMPLICIT DOUBLE FRECISION‘SA-HeF-2) 
ASIN = TATAN(X/DSQRT(1-x48%)) 

PRINTS», ° X =%s 2e * ASIN(X? 2°, ASIN 
Ke TURN 

END 


ro 


ee ee ee ee ee ee ee 


FUNCTION ACUS¢X) 

IMPLICIT DOUBLE FRECISION(A-HeF-Z) 
ACOS =~ TATANCLSQRT(1-X8X) 7X) 

FRINT®, ° x» ="s Xe * ACQS(X) =*s ACOS 
RE TURN 

EN! 
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FUNCTION EXFCX) 

IMPLICIT DOUBLE FRECISION(A-Hef-2) 
EXf = [WEXFUX) 

FRINT®, * X ="r Ky ° EXF(X) =". EXP 
RE TURN 

ENG 
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FUNCTION SQRKIT(X) 


90143 
00144 
00145 
00146 
v014? 


LMFLICIT) DOUBLE FRECISION(A-HeF-2) 
OSORTi 4) 


SQR' = 


PRINTS 
ht TUKN 
EN 


e °° X =%e 
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Q0149 
wUloe 
QuLot 
00152 
Oot 3 
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9Gb9¥ 
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JO160 
utist 


JUL. 
IQTAL 
00164 
0016S 


[ 


wD | 


) 


Me ° SORT, 


FUNCTION Alenia 


TPPOICT) DOUBLE FRECISIUN\A MeF-2) 
{iATaw 23 


“a AN - 
PRINTS 
he TURN 
tnt 


» * x 
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OOS 2 bBo CEOS T ARQ BE TUR” 
an Ss rat aet 2 LATAN ES AP a O57) ad eC COSNT/COST) 
‘Serge -{7.*N 
a EM. 
PLE S VERSIIN OO.) 


Ae ase 
VO262 FUNCTION FASINT CTHETAD 
09263 IMROICTY TOURLE PRECISTON(A-HrP-7} 


00264 EXTERNAL ASIN 


99265 COMAGN Nr Ab 

00246 FOST = DCOSC TRE TA: 
OOL4 ’ COSNT = OCOSCONSTHE TA) 
90268 FASINT . COSNI 


Qn &9 [fF <COST.EQ. 9) RETURN 
69279 FASINT - CASTNCARSCOST) AKI &CCOSNT/COST) 
R~20 


AD-Al116 070 GEORGIA INST OF TECH ATLANTA ENGINEERING EXPERIMENT =-ETC F/G 9/2 
INVESTIGATE CAPABILITY OF ADA HIGHER ORDER PROGRAMMING LANGUAGE-ETC (U) 
MAR 82 Ut J GALLAHER pry eerne gee 


UNCLASSIFIED RADC=-TR-82=-46 


2 ar 2 
AA 
1a, Pe 


00271 
00272 


00273 
00274 
00275 
00276 
00277 
00278 
00279 
00280 
00281 
00282 


00283 
00284 
00285 
00286 
00287 
00288 
00289 
00290 


00291 
00292 
00293 
00294 
00295 
00296 
00297 
00298 
00299 
00300 
00301 
00302 


00303 
00304 
00305 
00306 
00307 
00308 
00309 
00310 
00311 


RETURN 
END 
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FUNCTION FACOST( THETA) 

IMPLICIT DOUBLE PRECISION(A-HrP-2) 
EXTERNAL ACOS 

COMMON No AK 

COST = DCOS(THETA) 

COSNT = DCOS(NSTHETA) 

FACOST = ACOS(AK&COST) SCOSNT 
RETURN 

END 
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FUNCTION RASIN(X) 

IMPLICIT DOUBLE PRECISION(A-H»P-Z) 
EXTERNAL ASIN 

RASIN = 1/ASIN(X) 

PRINTS» ° X =", Xe ° RASIN(X) =" RASIN 
RETURN 

END 
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FUNCTION FRASINT (THETA) 

IMPLICIT DOUBLE PRECISION(A-HrP-Z) 
EXTERNAL RASIN 

COMMON Ns AK 

COST = DCOS(THETA) 

COSNT = DCOS(NSTHETA) 

FRASINT = COSNT 

IF (COST.EQ.0) RETURN 

FRASINT = COST&RASIN(AKSCOST ) SAKSCOSNT 
RETURN 

END 
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FUNCTION FALNT (THETA) 

IMPLICIT DOUBLE PRECISION(A-H»P~Z) 
EXTERNAL ALN 

COMMON Ny AK 

COST = DCOS( THETA) 

COSNT = DCOS(NSTHETA) 

FALNT = COSNT 


IF (COST.EQ.0) RETURN 
B-21 


00312 
00313 
00314 


QO322 
00323 
00324 
00325 
00324 
00327 
00338 
00329 
O034% 


NOS51 
00332 
0035223 
00334 
OO325 
99335 
00337 
00336 
003239 


90346 
O74) 
gagar 
Q0545 
G0344 
00345 
90546 
0034?) 
OO34E 


FALNT = (ALN(AKSCOST)/AK) &(COSNT/COST) 
RETURN 
END 
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FUNCTION ALN(Z) 

IMPLICIT DOUBLE FRECISION(A-HrF-2) 
ALN = PLOGC(Z + 19/02 - 20072 
FRINT&, © Z =*, Zs * ALN(Z? =*) ALN 
RETURN 

END 
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FUNCTION FLOGIC IHL TA) 

IMPLICIT DGQUELE FRECLS {ON (A-HeP~-Z) 
COMMON Ns AKs fy AL 

COST = PCOS(THETA: 

COSN) = LOCUS (NXTHE TA: 


FLOGT = DLOGC(AR + AL t CAR - ALD &COST); 27 kCOSNT 


RE TURN 
FN]: 
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FUNCTION FSQRTT (THETA) 

IMPLICIT DGUELE FRECTCION (A-HrF-Z) 
COMMON He Ay Ty AL 

POST = DeNSt THETA) 

COSNF «Tet eh HE TAD 


FSARTT = US@RTCCAR + 4. & CAR - ALD &COST)/2)&COSNT 


RETURN 
Eiti 
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FUNCTIUIN FEXPTCTHE TAD 

IMFLICL) SUUBLE URECISION (A-HeF-Z) 
COMMON Ny ARy J» Al 

COST = PCOS: THETA: 

COSNT = HCOS¢NATHE TA! 


FEXFT = DEXP( CAR + AL 4 CAR - ALD RCOST)/2)8COSNT 


RETURN 
ENII 
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00349 C 


00350 
00351 
00352 
00353 
00354 
00355 
00356 
00357 
00358 


00359 
00240 
00361 
00342 
003482 
00364 
60365 
GO5Sé 


00367 
00368 
00369 
00370 
00371 
00372 
00373 
00374 
00375 


00376 
00377 
00378 
00379 
00380 
00381 
00382 
00383 
00384 
00385 
00386 
00387 


00388 


FUNCTION FAASNT( THETA) 

IMPLICIT DOUBLE FRECISION(A-H»P-2) 

EXTERNAL ASIN 

COMMON Ne AKs Is AL 

COST = DCOS( THETA) 

COSNT = DCOS'N&THETA) 

FAASNT = ASING (AP + AL + (AR - AL) SCOST)/2)8COSNT 
RE TURN 

EN 


(FLECS VERSION 22.91) 


Eee Seer ewt ewe ee 2 ew ewe Se TF ert wwe eT BS eS eS OSB SO esr ew ST Oe se 


FUNCTION EX1X(X) 
IMPLICIT DOUBLE FRECISION (A-HrP-Z) 
EX1X = 1 
IF (X.NE.O) EXiX = (DEXP(X) - 1)/X 

C PRINT Rs ° X ="s ky ° EXIX(X) ="s EXUX 
RE TURN 
ENL 
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FUNCTION FEXLXT: THETA) 

IMPLICIT DOUBLE PRECIS1UN (A-HeP-Z) 

COMMON No AKys IT» AL 

CUST = DCOS< THETA) 

COSNT = DCOS(NXTHETA? 

FEXIXT = EX1X(CAK + AL + CAK - ALD SCOST)/2)SCOSNT 
RETURN 

END 
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FUNCTION SQAC(Y) 
C NOTE . .« THIS FUNCTION COMFUTES 1/2 OF ACOS(1-Y)aK2,/'Y 
IMPLICIT DOUBLE PRECISTION(A-HsP-Z) 
C EXTERNAL ACOS 
SQAC = t 
IF (Y.EQ.9) RETURN 
C SQAC = ACOS(L - Y)RR2/(QBY) 
SQAC = DATANCDSQRT(YR(2 — YoI/CL - ¥))882/(028Y) 
C PRINT®s “Y = *s Ys * SQACCY) = °s SQAC 
RETURN 
END 
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00390 
00391 
00392 
00393 
00394 
00395 
00396 


00397 
003598 
00399 
00400 
00401 
00402 
00403 
00404 


00405 
00406 
00407 
90408 
00409 
00410 
00411 
vo4al? 


00413 
00414 
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00417 
00418 
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00420 
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00423 
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00426 
00427 
00428 
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FUNCTION FSQACT( THETA) 

IMPLICIT DOUBLE PRECISION(A-HeF'-Z) 

COMMON Ne AKe Ie AL 

COST = DCOS( THETA) 

COSNT = DCOS(NSTHETA) 

FSGACT = SQAC( (AK + AL + (AR - AL) BCOST)/2)SCOSNT 
RETURN 

END 
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FUNCTION SINH(X) 

IMPLICIT GQOUBLE FRECISION(A-HeF-Z) 
E = DEXP(X) 

SINH = (E - 1/E)/2 
PRINT®s ° X = "so Xs ° 
RETURN 

ENT 


SINH(X) = °s SINH 
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FUNCTION COSH(X: 

IMPLICIT DOUBLE FRECISION(A-HeP-Z) 
E = JiEXP< Xx) 

COSH = \E t+ LEDS 2 
PRINTKs *° K = "so he ° 
RE TURN 

ENL 


COSH\X) = °+ COSKH 
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FUNCTION FSINHT: THETA, 

IMPLICIT DOUBLE FRECISION(A-HeF-Z) 
EXTERNAL SINH 

COMMON Ne AK 

COs = NCOSi THETA) 

COSN? = DCIS(NSTHETA) 

FSINHT = LOSNT 

IF (COST.FRQ.0) NETIRN 

FSINHT = (SINH(AKELUS, ). AK) AC COSNT/COST) 
RETURN 

ENT 
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FUNCTION FCSHT THETA) 
(MPLICIT GOQUBLE FRECISTONCA-HeF-2Z)} 
EXTERNAL COSH 


COMFON Ny At B-24 


00430 
00431 
00432 
00433 
00434 


00445 
OU446 
00447 
00448 
00449 
09450 
0451 
09432 
004353 
00434 
00455 
00456 


YO4S? C 


094353 
VUASY 
Qo440 
JIVASI 
O45? 
O0453 
GO3S54 


004945 C 


VOIa4 
GO04S? 
00468 
9044° 
OG Pe 


COST = DCOS( THETA) 

COSNT = DCOS(NSTHETA) 

FCOSHT = COSH(AK&SCOST)&COSNT 
RETURN 

END 
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FUNCTION ASINH(X) 

IMPLICIT DOUBLE FPRECISION(A-H»eF-Z) 
EXTERNAL ATANA 

ASINH = 9.920 

IF (X .EQ@. 9) RETURN 

ASINH = ATANH:X/05QeT. 1 +& XKX)) 

PRIGTE: © X =",s Xp ° ASINH(X) ="s ASINH 
RETURN 

END 
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FUNCTION FASI&HIT THETA! 

IMPLICIT DOUBLE FRECISION(A Hort -Z) 
EXTERNAL ASINH 

CUMMON Ne AK 

COST = fCOSCTHETA: 

COSNT = GCOS ENKTHET A} 

FASANAT = COSNy 

{F «COST.EG 9, RETURH 

PASIGH? = CASLam eh eCOET: TAR ALCON  TOST) 
Ro TURAN 

EN Ii 
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FUNCTION ACOSiE <: 

IMFL i, © BOMKLE URRCTSION CAH S- 23 
EZTLENAL ATRNE 

ACOSH = ATANACTSQRTCRaK- 17a: 

PRINTKs “x 2%. 2 " AcuzHek) =" ACOSH 
RETURN 

ENII 
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EXTERNAL Altos 

COmawe tts ét 
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N04) 1 
0472 
00473 
Q9474 


V0479 


60476 
00477 
Y04/8 
10479 
OO480 


Og4el C 


00482 
00435 


OG4B4 + 


00485 
00434 
00487 
00488 
00489 
00470 
00491 
00492 
00493 
00494 
06495 


00496 
00497 
00498 
00499 
90500 
00501 
09502 


00503 
00504 
00505 
00506 
00507 
00508 
00509 
00510 
00511 


COSNT = DLUS (Ne SHETA? 
FACOSHT-= ACOSHCAKKCOST ) &COSNT 
RE TURN 

ENN 
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FUNCTION ATANH(X: 

IMPLICIT fOUBLE FREE CS CONCA~HeFH-Z) 
HTANH = 6.0106 

IF +X ,EQ. OF RETURN 

ATANH © BLUGSS14K)/¢1-xK)0/2 

PRINT¥s *© X :*5 Xe * ATANH(X) =°s ATANH 
RE TURN 

ENG 
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FUNCTION FATANHT( THETA) 

IMPLICIT DOUBLE FRECISION(A-HsF-Z) 
EXTERNAL ATANH 

COMMON Ne AK 

COST = DCOSCTHETA) 

COSNT = DCOQS¢NXTHETA) 

FATANHT = COSNT 

[IF (COST.EQ@.0) RETURN 

FATANHT = CATANH(AKXCOST)/AK) &(COSNT/COST) 
RETURN 

END 
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FUNCTION TANH(X) 

IMPLICIT DOUBLE PRECISION(A-HrP-Z) 
TANH = DTANH(X) 

PRINTS» ° X ="» Xe * TANH(X) =°s TANH 
RETURN 

END 
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FUNCTION FTANHT (THETA) 

IMPLICIT DOUBLE PRECISION(A-H»P-Z) 
EXTERNAL TANH 

COMMON Ny AK 

COST = BCOS( THETA) 

COSNT = DCOS(NXSTHETA) 

FTANHT = COSNT 


IF (COST.EQ@.0) RETURN 
B=-26 


00512 FTANHT = (TANHCAK&COST)/AK) &(COSNT/COST) 
00513 RETURN 
00514 END 


| 
! 
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EOI. O FILES. 1 RECS. 2321 WORDS. 
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99997 
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9999 / 


FUNCTION ROMBRG(A» By Fe OE) 
IMPLICIT DOUBLE FRECISION(A-HeP-Z) 
INTEGER FOUR(21)» 0 

DIMENSION Q(20)* Q0(2) 

EQUIVALENCE (Q0(2)98(1)) 

DATA 0/0/» LMAX/1/s FOURC1)/4/ 


Q(0) = 0 
H = (BR - A)/2 
IMAX = 1 


BO 99999 L=1720 
IF(.NOT.(L.»GE-LMAX)) GO TO 997998 
LMAX =L +1 
FOUR(L+1) = 4FOQUR(L) 
SUM = 0 
DO 99997 T=1*9IMAX 
= (2kI - IMAX - 1.020) /1MAX 
= UKCS - UEYV)/2 
Fo= 1,520k(1 - Vey) 
X=A +t (U + 1)KH 
SUM = SUM + F(X) ®UF 
CONTINUE 
SUM = SUM + IMAXKQ(L-1) 
IMAX = 2KIMAX 
Q¢Ly = SUM/IMAZ 
10 79994 M=islk. 
K=t-t 
QR = CFOURCHIKQCKE1S - QIK) /SCFOUR(M) - 13 
CONTINUE 
IF(.NOT.¢L.GT.1)7 GG TO 99995 
ROMERG = (B - AdKDID) 
QQ - fARSEQCO) - GED KCBR - Ad? 
[F (QTIQ,LE.0OE) RETURN 
UNE ENUE 
CONTINUE 
RE TURN 
ENL 
FUNCTION ROMBER CAs Be Fy OE) 
IMFLICTT DOUBLE FRECISION(A-HeF-Z) 
INTEGER FOURS 21)» O 
DIMENSION (204 QOC2) 
EQUIVALENCE (QG°62)sQC10) 
ATA O/0/9 LMAX7i*ss) FOURCT 747 
Ni) = FCA) + FCBY)/2 
H:- K-~- A 
lax = J 
Itt) 9999G La19 00 
TFC .NOT,(L.GE.L MAX)? GO TO 99998 
(MAX = #1. 
FOUR ‘L419 = 4kFOURCL) 
H: H/2 
SUM. 9 
DO) 99997 T=1+ {MAX 
SUM = Sum + F(A + HRC2KT - 1)) 
CUNT ENUE 
SUM - SUM + TmMAaxkQsl-1) 
[MAX = 2XIMAT 
Q¢L> = SUM/IMAX 
LO 99996 M=Lel 
K=l- 
ACK) = (FOURCM)SQCKEL) - ae - 1) 


V 
UY 
U 


60008 
00009 
00011 
00012 
60013 
00015 
99017 
00018 
00019 
00021 
06022 
00023 
00624 
09027 
00628 
00029 
00030 
00031 
00032 
00033 
90034 
00035 
00036 
00037 
00039 
00046 
00041 
00042 
00044 
00045 
00050 
00051 
00054 
00054 
00055 
00056 
00063 
000464 
00066 
0004? 
00068 
00070 
09072 
00073 
00074 
00076 
00077 
00073 
00079 
00082 
00083 
00084 
00084 
00084 
00085 
00084 
00087 
00089 
00090 
00091 


99996 CONTINUE 


99999 


ROMBER = (B - &)#Q(0) 

ODQ = DABS((Q(0) - Q(1))&(B - A)) 
IF (ODQ.LE.OE) RETURN 

CONTINUE 

RETURN 

END 

FUNCTION SIN(X) 

IMPLICIT DOUBLE FRECISION(A-H»P-Z) 

SIN = BSIN(X) 

RETURN 

EN 

FUNCTION TAN(X) 

IMPLICIT ROUBLE FRECISION(A-HsP-Z) 
TAN = DSIN(X)/DCOS(X) 

RETURN 

ENE 

FUNCTION ASIN(X) 

IMPLICIT ROUBLE PRECISION( A-H»P-Z) 
ASIN = DATAN(X/DSQRT(1-X%X)) 
RETURN 

ENL 

FUNCTION ACOS(X) 

IMFLICIT BOUBLE FRECISION(A-HsF-Z) 

ACOS = DATANCDSQRT(1-X*X)/X) 

RETURN 

ENT 

FUNCTION EXFC(X) 

IMPLICIT DOUBLE FRECISION(A-H»P-Z) 

EXF = DIEXP(X) 

RETURN 

ENE 

FUNCTION SQRT(X) 

IMFLICIT BOUBLE FRECISION(A-H»P-Z) 
SQRT = DSQRT(X) 

RETURN 

END 

FUNCTION ATAN(X) 

IMFLICIT DOUBLE PRECISION(A-H»F-Z) 

ATAN = TATAN(X) 

RETURN 

ENE 

FUNCTION ALQG(X) 

IMPLICIT DOUBLE FRECISION(A-H»P~Z) 
ALOG = DPLOG(X) 

RETURN 

END 

FUNCTION COS(X) 

IMPLICIT DOUBLE FRECISION( A-H»P-Z) 
COS = DCOS(Xx) 

RETURN 

END 

FUNCTION FCOST (THETA) 

IMPLICIT DOUBLE PRECISION (A-H,P-Z) 

COMMON Ne AK 

COST = DCOS( THETA) 


COSNT = DCOS(NKTHETA) 

FCOST = DCOS(AKXCOST) XCOSNT 
RETURN 

END 


FUNCTION FSINT(THETA) B-29 


00092 
00094 
00099 
00100 
00102 
00103 
00104 
00107 
00108 
00109 
00111 
00112 
00114 
00115 
09116 
00118 
OO119 
00121 
0012? 
O0123 
06125 
00126 
006128 
00129 
001390 
00132 
00133 


00163 
00164 
00165 
0016? 
00148 
00170 
00171 
00172 
00173 
00174 
00175 
00176 
00177 
00179 


99999 


99597 


79996 
9999? 
99998 
99999 


IMPLICIT DOUBLE PRECISION(A-H»F-Z) 
COMMON Ny AK 

COST = fCOS<« THETA) 

COSNT = [ICOS(NXTHETA) 

FSINT = COSNT 

IF <«COST.EQ.0O) RETURN 

FSINT = (CUSINCAK&COST) /AN) *(COSNT/COST) 
RETURN 

END 

FUNCTION FOLY(X» Ar N) 

IMPLICIT DOUBLE FRECISION (A-HeF-Z) 
TIMENSION AC25) 

PROD = ACN) 

[10 99999 J=1l9N 

PROG = PRODUXX + ACN-J) 

CONTINUE 

FOL, = RROD 

RE TURN 

ENS) 

SUBROUTINE ACOEF CAs By ICs AKy N) 
IMPLICIT DOURLE FRECISTON( A-HeF-7) 
CIMENSION BelSbts ACLSts WES oS: 
lu Y999y J=1leN 

SUM = vu 

LO 99993 h=JeNn 

Sum =SUM ¢ BCKIKICKR: J) 


; CUNT INUE 


nid? = SUM/AKRK 

CUNT ENUE 

bg 

hK = 0 

SUM = BIR GRICtK® J) 

Pu 9Se97 poten 

wilh = sum + Bats xIConwe OD 

ONY Tale 

Aide = SUM 

RETURN 

ENT 

SUBROUTINE TCH ICs a 

LTMENSIUN [CG OSs CS: 

INTEGER U 

fiATé ay J, 

FOtie: We 

IC(tt-le DQ) - y 

HD Poses RAL -N 

Toths Qs = -EC«K-29 OQ 

LOGh ee Dt SMe 

\CiQ-t.s 2) + | 

LU 79995 J=1eh 

TFe. POT. MG Jthe 2) EQ.39) GO TO 29995 

Thite Pot a 

GA i 89O97 

TOKs Joe UHLE Oe Le let) - Winete ot 

CONT LNUE 

CONTINUE 

PONT INUE 

RETURN 

ENT 

FUNCTION FTANT. THETA? 

UMPL ICE E LOU LE PRECISLUN(A-HeP-2) 
CTERNAL B-30 


00180 
00181 
00182 
00183 
00184 
00185 
00184 
00167 
00183 
00191 
06192 
69193 
001974 
OOl7S 
O97 
00195 
G919S 
OU19? 
O919% 
Qunog 
Cy ad 
Ouced 
Yyl)4 
OO205 
OO 
we CUS 


bot 


oe 


~ 


(Gr 


~ 
‘J 
ache ? 
a, 


ed 


COMMUN the AK 


bast = bschHe la? 
PiSNE > LCOUSONS THETA 
kiatl LUSH 


LF FUN cleus RETURN 

FANT = CP ANCARACOST) ZAR a CLUSNT/COST) 
RE TUR! 

tat 

Fite PTIN FATANT OTHE TA) 

TMFLIL EY DOUBLE FRECISLONCA-H»F-Z) 
Cmetiks Ne AK 


most o- feast YHE TA) 

Citsht - 00S \deTHET A} 

FATAN) = CUSNE 

It CET EW.) RETURN 

FATANT = (DAIANCARACUST AK RICOSNT/COST) 
ht TURN 

ENT 


FUNCTION FASINT< THETA) 

PMFLIGL) DOUBLE FRECLSTON(A-HeF-Z) 
EXTERNAL ASIN 

LUMMUN Ne» AN 

COST = fCUSCTHETA} 

COSNT = DCOSCNETHETA) 

FASINT = COSNT 

TF «(eOST.EQ.0% RETURN 

FASTNT = (ASIN(AKECOST) /AK) K(COSNT/COST > 
RETURN 

ENT 

FUNCTION FACOST( THETA) 

TMPLICIT TNOURLE FPRECISION(CA-HeF-Z) 
EXTERNAL ACOS 

COMMON Ne AK 

COST = DCOS(THETA) 

COSNT = DCOS(NXTHETA) 

FACOST = ACOS(CARKKCOST) XCOSNT 
RETURN 

END 

FUNCTION RASIN(X: 

IMPLICIT DOUBLE FRECISION(A-HrF-Z) 
EXTERNAL ASIN 

RASIN = 1/ASIN(X) 

RETURN 

END 

FUNCTION FRASINT (THETA) 

IMPLICIT DOUBLE FRECISION(A-HsP-Z) 
EXTERNAL RASIN 

COMMON Ny AK 

COST = DCOS( THETA) 

COSNT = DCOS(NXTHETA) 

FRASINT = COSNT 

IF (COST.EQ.0) RETURN 

FRASINT = COSTKRASIN(AKXCOST )*AKKCOSNT 
RETURN 

END 

FUNCTION FALNT (THETA) 

EMPLICIT fOUBLE FRECISION(A-H»F-Z) 
EXTERNAL ALN | 
COMMON Ne AK 

COST = DCOS(THETA) 


COSNT = DCOS(NETHETA) B-31 


00242 
00243 
00244 
00245 
00246 
00247 
00248 
00249 
00251 
00252 
00253 
00254 
00255 
00254 
00257 
00258 
00259 
00260 
00262 
00263 
00264 
00265 
00266 
00267 
00268 
60269 
00270 
00271 
00272 
00274 
00275 
00276 
00277 
00279 
00279 
00280 
00281 
00282 
00284 
00285 
00286 
00287 
00289 
00290 
00292 
00293 
00294 
00295 
00296 
00297 
00298 
00299 
00300 
00301 
00302 
00304 
00305 
00306 
00307 
00308 
00309 


FALNT = COSNT 

IF (COST.EQ.0) RETURN 

FALNT = (ALN(AKACOST)/AK)&(COSNT/COST) 
RETURN 

END 

FUNCTION ALN(Z) 

IMPLICIT GOUBLE FRECISION(A-HsF -2) 
ALN = [ILOG((7 + 13/01 - Z))/2 
RETURN 

END 

FUNCTION FLOGT( THETA) 

IMPLICIT MOUBLE FRECISION (A-HeF-Z) 
COMMON Ny AKy Ivy AL 

COST = [COS( THETA) 


COSNT = DICOS(NXTHETA) 

FLOGT = BLOG (AK + AL + (AK - AL) &COST)/2)*COSNT 
RETURN 

ENT 


FUNCTION ¢SQRTT( THETA) 

IMPLICIT DOUBLE PRECISION (A-H»F-Z) 
COMMON Ns AKy Is AL 

COST = DCOS( THETA) 

COSNT = DCOS(NXTHETA) 

FSQRTT = BPSQRTCCAK + AL + (AK - AL) *COST)/2)xCOSNT 
RETURN 

ENT 

FUNCTION FEXFT¢( THETA? 

IMFLICTT DOUBLE FRECISION (A-H»yF-Z) 
COMMON Ny AKy Ie AL 

COST = fICOS( THETA) 

COSNT = DCOS<(NXTHETA) 

FEXFT = DEXFCCAK + AL + (AK - AL) *XCOST)/2)&COSNT 
RETURN 

END 

FUNCTION FAASNTI THETA) 

IMPLICIT (iGUBLE FRECISION(A-HoF-Z) 
EXTERNAL ASIN 

COMMON Ny AKy Is AL 

COST = DCOS( THETA? 

CQOSNT = f1COS(NKTHETA) 

FAASNT = ASIN(CCAK + AL + (AK - AL) *COST)/2)xCOSNT 
RETURN 

ENT! 

FUNCTION EX1*%¢(x) 

TMFLICIT fOUBLE FRECISION (A-HsF-Z) 
Exia = | 

IF (X,NE,O; EX1xX = (DEXFOCX) - 1)/% 
RETURN 

ENT! 

FUNCTION FEX1IXT¢ THETA) 

IMFLICI} [OQUBLE PRECISION (A-H»yF-Z) 
COMMON Nes AKs Is AL 

COST = DCOS( THETA) 

COSNT = [ICOS(NRTHETA) 

FEXiXT = EXIx( (AK + AL + (AK - AL) ECOST)/2)*COSNT 
RETURN 

ENLi 

FUNCTION SQAC(Y) : 
IMPLICIT DOUBLE FRECISION(CA-HsF-Z) 
SQAC = 1 


IF (1.E0.0) RETURN ae 


00310 
00311 
00312 
00313 
00314 
00316 
00317 
00318 
00320 
00321 
00323 
90324 
00325 
00326 
00327 
00328 
00329 
00330 
00332 
00333 
0035354 
00335 
00336 
00337 
00338 
00339 
00341 
00342 
00343 
00344 
003545 
00346 
00347 
00348 
00350 
00351 
00352 
003353 
00354 
00355 
00354 
00357 
00358 
00360 
00361 
00362 
00363 
00365 
00366 
00348 
00369 
00370 
00371 
00372 
00373 
00374 
00375 
00377 
00379 
00381 
00382 


SQAC = DATAN(IISQRT¢(YR(2 - Y))/C1 ~- Y))&R2/(28Y) 


RETURN 

END 

FUNCTION FSQACT( THETA) 

IMPLICIT fQUBLE FRECISION(A-HeF-Z) 
COMMON Nr AKy Ie AL 

COST = DCOS( THETA) 

COSNT = QUOS¢NEXTHETA) 


FSQACT = SQAC(CAK + AL + (AK - AL) COST) /2)xCOSNT 


RETURN 

ENE 

FUNCTION SINH(X) 

IMPLICIT fGUBLE FRECISION(CA-HeF-Z) 
£ = PEXFCX) 

SINH = «—£ ~- L/Ed 72 

RETURN 

ENT 

FUNCTION COSH: x) 

iMFLICIT DOUBLE FRECISION(A-HeF-Z) 
E = TEXFCX) © 

COSH = (E + 1vEs 72 

RETURN 

ENLt 

FUNCTION FSINHTi THETA) 

IMPLICIT {CUBLE FRECISTON(A-HeF-Z) 
EXTERNAL = fNH 

COMMUN Ne Ah, 

CUST = RCOS:THETA: 

COSNT = LCOS CN#THETA) 

FSINHT = COSNT 

IF (COST.EQ.0: RETURN 

FSINH! = (CSLNHCARKCUST) SAR HCLOUSNT/COSY 
RETURN 

ENT 

FUNCTION FCOSHT (THETA) 

IMPLICIT DOQUELE FRECISION(A-HsF-2Z)> 
EXTERNAL COSH 

COMMON Ne AK 

COST = COS (THETA? 

COSNT = LCOS (NKTHETA) 

FCOSHT = COSH(AKKCOST) &COSNT 
RETURN 

END 

FUNCTION ASTNH(X) 

IMFEJCIT TOUBLE FRECISION(A-H»F-Z) 
EXTERNAL ATANH 

ASINH = 0.0[0 

IF (X .EQ, OF RETURN 

ASINH = ATANH(X/ISQRT(1 + XEX)) 
RETURN 

ENT 

FUNCTION FASINHT*‘ THETA) 

IMPLICIT DOUBLE FRECISION(CA-HsF-Z) 
EXTERNAL ASINH 

COMMON Ne AK 

FOST = MCOS( THETA) 

COSNT = fICOS(NXTHETA) 

FASINHT = COSNT 

IF (COST.EQ.0) RETURN 


FASINHT = (ASINH(AKSCOST) /AK) &(COSNT/COST) 


RETURN 
B-33 
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MISSION 
of 
Rome Air Development Center 


RADC plans and executes research, development, test and 
selected acquisition programs in support of Command, Contnrok 
Communications and Intelligence (C31) activities. Technical 
and engineering support within areas of technical competence 
4& provided to ESP Program Offices (P08) and other ESP 
eLements. The principal technical mission areas are 
communications, electromagnetic guidance and controk, sur- 
veiklance of ground and aerospace objects, AntellLigence data 
colfection and handling, information system technology, 
Aonospheric propagation, solid state ACLENCLSA, microwave 
physics and electhonic reliability, marintainability and 
compatibility. 
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